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Abstract 

The four LEP collaborations, ALEPH, DELPHI, L3 and OPAL, have searched for the neutral 
Higgs bosons which are predicted by the Minimal Supersymmetric Standard Model (MSSM). 
The data of the four collaborations are statistically combined and examined for their consistency 
with the background hypothesis and with a possible Higgs boson signal. The combined LEP 
data show no significant excess of events which would indicate the production of Higgs bosons. 
The search results are used to set upper bounds on the cross-sections of various Higgs-like event 
topologies. The results are interpreted within the MSSM in a number of "benchmark" models, 
including CP-conserving and CP-violating scenarios. These interpretations lead in all cases to 
large exclusions in the MSSM parameter space. Absolute limits are set on the parameter tan f3 
and, in some scenarios, on the masses of neutral Higgs bosons. 
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1 Introduction 



One of the outstanding questions in particle physics is that of electroweak symmetry breaking 
and the origin of mass. The leading candidate for an answer is the Higgs mechanism [1] whereby 
fundamental scalar Higgs fields acquire nonzero vacuum expectation values and spontaneously 
break the electroweak symmetry. Gauge bosons and fermions obtain their masses by interacting 
with the resulting vacuum Higgs fields. Associated with this description is the existence of 
massive scalar particles, the Higgs bosons. 

The Standard Model [2] requires one complex Higgs field doublet and predicts a single 
neutral Higgs boson of unknown mass. After extensive searches at LEP, a lower bound of 
114.4 GeV/c 2 has been established for the mass of the Standard Model Higgs boson, at the 
95% confidence level (CL) [3]. 

Supersymmetric (SUSY) [4] extensions of the Standard Model are of interest since they pro- 
vide a consistent framework for the unification of the gauge interactions at a high energy scale 
and for the stability of the electroweak scale. Moreover, their predictions are compatible with 
existing high-precision data [5]. The Minimal Supersymmetric Standard Model (MSSM) (re- 
viewed, e.g., in [6]) is the SUSY extension with minimal new particle content. It requires two 
Higgs field doublets and predicts the existence of three neutral and two charged Higgs bosons. 
The lightest of the neutral Higgs bosons is predicted to have a mass less than about 140 GeV/c 2 
including radiative corrections [7] . This prediction provided a strong motivation for the searches 
at LEP energies. 

Most of the experimental investigations carried out in the past at LEP and elsewhere were 
interpreted in MSSM scenarios where CP conservation in the Higgs sector was assumed. In 
such scenarios the neutral Higgs bosons are CP eigenstates. However, CP violation in the 
Higgs sector cannot be a priori excluded [8]. Scenarios with CP violation are theoretically 
appealing since they provide one of the ingredients needed to explain the observed cosmic 
matter-antimatter asymmetry. The observed size of CP violation in B and K meson systems is 
not sufficient to drive this asymmetry. In the MSSM, however, substantial CP violation can be 
induced by complex phases in the soft SUSY-breaking sector, through radiative corrections, es- 
pecially from third-generation scalar quarks [9] . In such scenarios the three neutral Higgs mass 
eigenstates are mixtures of CP-even and CP-odd fields, with production and decay properties 
different from those in the CP-conserving scenarios. Hence, the experimental exclusions pub- 
lished so far for the CP-conserving MSSM scenarios may be weakened by CP-violating effects. 
There is currently one publication on searches interpreted in CP- violating scenarios [10]. 

In this paper we describe the results of a statistical combination based on the searches of the 
four LEP collaborations [10-13], which was carried out by the LEP Working Group for Higgs 
Boson Searches. These searches include all LEP2 data up to the highest energy, 209 GeV; in 
the case of Refs. [10, 12] they also include the LEP1 data collected at energies in the vicinity 
of 91 GeV (the Z boson resonance). The combined LEP data show no significant signal for 
Higgs boson production. The search results are used to set upper bounds on topological cross- 
sections for a number of Higgs-like final states. Furthermore, they are interpreted in a set of 
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representative MSSM "benchmark" models, with and without CP-violating effects in the Higgs 
sector. 



2 The MSSM framework 

The LEP searches and their statistical combination presented in this paper are interpreted in 
a constrained MSSM model. At tree level, two parameters are sufficient (besides the known 
parameters of the Standard Model fermion and gauge sectors) to fully describe the Higgs sector. 
A convenient choice is one Higgs boson mass {mx is chosen in CP-conserving scenarios and mn± 
in CP- violating scenarios), and the ratio tan/5 = v^/vi of the vacuum expectation values of 
the two Higgs fields (v 2 and v\ refer to the fields which couple to the up- and down-type 
fermions). Additional parameters, MgusY; M 2 , fi, A and mg, enter at the level of radiative 
corrections. Msusy is a soft SUSY-breaking mass parameter and represents a common mass 
for all scalar fermions (sfermions) at the electroweak scale. Similarly, M 2 represents a common 
SU(2) gaugino mass at the electroweak scale. The "Higgs mass parameter" \i is the strength of 
the supersymmetric Higgs mixing; A = A t = is a common trilinear Higgs-squark coupling 
at the electroweak scale and mg the gluino mass. Three of these parameters define the stop and 
sbottom mixing parameters X t = A — /i cot (3 and — A — \l tan f3. In CP- violating scenarios, 
the complex phases related to A and mg, arg(A) and arg(mg), are supplementary parameters. 
In addition to all these MSSM parameters, the top quark mass also has a strong impact on 
the predictions through radiative corrections. In this paper, four fixed values are used in the 
calculations: m t = 169.3, 174.3, 179.3 and 183.0 GeV/c 2 . For the purposes of illustration, m t 
= 174.3 GeV/c 2 is used in producing the figures (unless explicitly specified otherwise), which 
is a previous world-average value [14] and which is within the current experimental range of 
172.7±2.9 GeV/c 2 [15]. The influence of the top quark mass on the exclusion limits is discussed 
in Sections 5 and 6 along with the other results. 

The combined LEP data are compared to the predictions of a number of MSSM "bench- 
mark" models [16]. Within each of these models, the two tree-level parameters, tan/3 and 
tua (in the CP-conserving scenarios) or mn± (in the CP-violating scenarios) are scanned while 
the other parameters are set to fixed values. Each scan point thus represents a specific MSSM 
model. The ranges of the scanned parameters and the values of the fixed parameters are listed 
in Tabled for the main scenarios studied. The first five models represent the main benchmarks 
for CP-conserving scenarios while the last model, labelled CPX, is a benchmark model for CP- 
violating scenarios. Some variants of these benchmark scenarios, which are also investigated, 
are presented in the text below. 

The scan range of tan f3 is limited by the following considerations. For values of tan /3 below 
the indicated lower bounds, the calculations of the observables in the Higgs sector (masses, 
cross-sections and decay branching ratios) become uncertain; for values above the upper bounds, 
the decay width of the Higgs bosons may become larger than the experimental mass resolution 
(typically a few GeV/c 2 ) and the modelling of the kinematic distributions of the signal becomes 
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Benchmark parameters 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






no-mixing 


large-fx 


gluophobic 


small-a e ff 


CPX 


Parameters varied in the scan 


tan (3 


0.4-40 


0.4-40 


0.7-50 


0.4-40 


0.4-40 


0.6-40 


m A (GeV/c 2 ) 


0.1-1000 


0.1-1000 


0.1-400 


0.1-1000 


0.1-1000 




m H ± (GeV/c 2 ) 












4-1000 


Fixed parameters 


Msusy (GeV) 


1000 


1000 


400 


350 


800 


500 


M 2 (GeV) 


200 


200 


400 


300 


500 


200 


H (GeV) 


-200 


-200 


1000 


300 


2000 


2000 


TO| (GeV/c 2 ) 


800 


800 


200 


500 


500 


1000 


X t (GeV) 


2 Msusy 





-300 


-750 


-1100 


A — n cot (3 


A (GeV) 


X t +/j,cot P 


X t +/J,cot P 


X t +/j,cot P 


Vt+^cot (5 


X t +/j,cot (5 


1000 


arg(^4)=arg(m g ) 












90° 



Table 1: Parameters of the main benchmark scenarios investigated in this paper. The values 
of tan/5 and the mass parameters m A (in the CP-conserving scenarios) or mn± (in the CP- 
violating scenarios) are scanned within the indicated ranges. For the definitions of A and 
X t , the Feynman-diagrammatic on-shell renormalisation scheme is used in the CP-conserving 
scenarios and the MS renormalisation scheme in the CP-violating scenarios. 

inaccurate 2 . The scan range of m A is limited in most cases to less than 1000 GeV/c 2 ; at higher 
values the Higgs phenomenology is insensitive to the choice of m A . 

For a given scan point, the observables in the Higgs sector are calculated using two theo- 
retical approaches, both including one- and two-loop corrections. The FeynHiggs2.0 code [17] 
is based on a Feynman-diagrammatic approach and uses the on-shell renormalization scheme. 
The SUBHP0LE calculation and its CP-violating variant CPH [18] are based on a renormalization- 
group improved effective potential calculation [19] and use the MS scheme 3 . 

In the CP-conserving case, the FeynHiggs calculation is retained for the presentation of the 
results since it yields slightly more conservative results (the theoretically allowed parameter 
space is wider) than SUBHP0LE does. Also, FeynHiggs is preferred on theoretical grounds since 
its radiative corrections are more detailed than those of SUBHP0LE. 

In the CP-violating case, neither of the two calculations is preferred on theoretical grounds. 
While FeynHiggs contains more advanced one-loop corrections, the CPH code has a more precise 
phase dependence at the two-loop level. We opted therefore for a solution where, in each scan 
point, the CPH and FeynHiggs calculations are compared and the calculation yielding the weaker 

2 The DELPHI Collaboration included the variation of the Higgs boson decay width with tan/3 in their 
simulation for tan/3 between 30 and 50. With increasing tan/3, DELPHI observed an increase of the mass 
resolutions and hence a loss in the signal detection efficiencies; but this was compensated by the increase of the 
cross-sections, such that DELPHI found no significant drop in the overall sensitivity. 

3 New developments in this approach are implemented in the code CPsuperH [20]. 
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exclusion (more conservative) is retained. However, we also discuss in Section 6 the effect of 
using separately either one or the other of the two calculations. Rather large discrepancies 
between the two codes are found in calculating the partial width for the Higgs boson cascade 
decay T(H2 — ► H.1H1) (Hi and H2 are the lightest and the second-lightest neutral MSSM Higgs 
bosons). Aiming at conservative exclusion limits, therefore, the CPH formula for this decay was 
also used within the FeynHiggs code. 

All codes are implemented in a modified version of the HZHA program package [21], which 
takes into account initial-state radiation and the interference between identical final states from 
Higgsstrahlung and boson fusion processes. 

2.1 CP-conserving scenarios 

Assuming CP conservation, the spectrum of MSSM Higgs bosons consists of two CP-even 
neutral scalars, h and H (h is defined to be the lighter of the two), one CP-odd neutral scalar, 
A, and one pair of charged Higgs bosons, H 1 * 1 . The following ordering of masses is valid at tree 
level: m h <(M z , mjCjKm-a and m w ±<m H ±. This ordering may be substantially modified by 
radiative corrections [7] where the largest contribution arises from the incomplete cancellation 
between top and scalar top (stop) loops. The corrections affect mainly the neutral Higgs boson 
masses and decay branching ratios. 

In e + e~ collisions at LEP energies, the main production processes of h, H and A are the 
Higgsstrahlung processes e + e~— > hZ and HZ and the pair production processes e + e~— > hA 
and HA (in most of the MSSM parameter space only the hZ and hA processes are possible 
by kinematics). The fusion processes e + e~— > (WW— > h)u c i7 c and e + e~— * (ZZ^ h)e + e~ play a 
marginal role at LEP energies but they are also taken into account in the derivation of the 
results. 

The cross-sections for Higgsstrahlung and pair production can be expressed in terms of the 
Standard Model Higgs boson production cross-section cr^. The following expressions hold for 
the processes involving the lightest scalar boson h: 

a hZ = sin 2 ^ - a) 4z (1) 

a hA = cos 2 ((3-a)\a™. (2) 

Here a is the mixing angle which diagonalises the CP-even Higgs mass matrix (at lowest order 
it can be expressed in terms of m^, M z and tan/?) and A is a kinematic factor: 

A = Ai?/[A^ 2 (12J^/ a + Azh)] (3) 

with 

Ay = [1 - (rrii + rrij) 2 /s][\ - (m* - m^/s], (4) 

where s is the square of the centre-of-mass energy. The cross-sections for the processes in- 
volving the heavy scalar boson H are obtained by interchanging the MSSM suppression factors 
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sin 2 (/3 — a) and cos 2 (/3 — a) in Eqs. 1 and 2 and replacing the index h by H in Eqs. 1, 2 and 3. 
The Higgsstrahlung and pair production cross-sections are complementary, as seen from Eqs. 1 
and 2. At LEP energies, the process e + e~— > hZ is typically more abundant at small tan/5 and 
e + e~— > hA at large tan/3, but the latter process can be suppressed also by the kinematic factor 
A. 

The following decay features are relevant to the neutral MSSM Higgs bosons. The h boson 
decays mainly to fermion pairs, with only a small fraction of WW* and ZZ* decays, since its 
mass is below the threshold of the on-shell processes h— > WW and h^ ZZ. However, for 
particular choices of the parameters, the fermionic final states may be strongly suppressed. 
The A boson also decays predominantly to fermion pairs, independently of its mass, since its 
coupling to vector bosons is zero at leading order. For tan/3>l, decays of h and A to bb 
and r + r~ pairs are preferred while the decays to cc and gluon pairs are suppressed. Decays 
to cc may become important for tan/3<l. The decay h— > AA may be dominant if allowed 
by kinematics [22]. Higgs boson decays into SUSY particles, such as sfermions, charginos or 
invisible neutralinos, are suppressed due to the high values of the SUSY-breaking scale Msusy 
which have been chosen. 

In the following we describe the CP-conserving benchmark scenarios [16] which are examined 
in this paper. The corresponding parameters are listed in Table ^ 

2.1.1 The m^-max scenario 

In the m^-max scenario the stop mixing parameter is set to a large value, X t = 2Msusy- This 
model is designed to maximise the theoretical upper bound on for a given tan j3 and fixed m t 
and Msusy (uncertainties due to unknown higher-order corrections are ignored). This model 
thus provides the largest parameter space in the direction and conservative exclusion limits 
for tan j3. 

We also examine a variant of this scenario where the sign of \i is changed to positive, since 
this is favoured by presently available results on (g — 2)^ [23,24]. This variant is labelled 
m^-max (a) below. Furthermore, we examine the case where, besides changing the sign of 
fi to positive, the sign of the mixing parameter X t is changed to negative. This choice of 
parameters gives better agreement with measurements of the branching ratios and of the CP- 
and isospin-asymmetries for the process b— > [16,25]. This variant is labelled m^-max (b) 
below. 

2.1.2 The no-mixing scenario 

In the no-mixing scenario the stop mixing parameter X t is set to zero, giving rise to a relatively 
restricted MSSM parameter space. We also examine a variant of this scenario where the sign of 
\i is changed to positive, for a better agreement with recent measurements of (g — 2)^ [23,24], 
and Msusy is raised to 2 TeV in order to enlarge the parameter space of the standard no-mixing 
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scenario [16]. This variant is labelled no-mixing (a) below. In this case, tan/3 is scanned only 
from 0.7 upward due to numerical instabilities at lower values in the diagonalisation of the mass 
matrix. 

2.1.3 Special scenarios 

Some scenarios were designed to illustrate choices of the MSSM parameters for which the 
detection of Higgs bosons at LEP, at the Tevatron and at the LHC is expected to be difficult 
a priori due to the suppression of some main discovery channels [16]. 

• The large- fi scenario is constructed in such a way that, while the h boson is accessible by 
kinematics at LEP for all scan points, the decay h^ bb, on which most of the searches 
at LEP and at the Tevatron are based, is typically strongly suppressed. For many of the 
scan points the decay h— > t + t~ is also suppressed, such that the dominant decay modes 
are h— > cc, gg and WW*. The detection of Higgs bosons thus relies mainly on flavour- and 
decay-mode-independent searches. Moreover, for some of the scan points, the e + e~— > hZ 
process is suppressed altogether by a small value of sin 2 (/9 — a). In such cases, however, 
the heavy neutral scalar H is within reach (m H < 111 GeV/c 2 ) and the cross-section for 
e + e~-^ HZ, proportional to cos 2 (/3 — a), is large; the search may thus proceed via the 
heavy Higgs boson H. 

• The gluophobic scenario is constructed in such a way that the Higgs boson coupling to 
gluons is suppressed due to a cancellation between the top and the stop loops at the hgg 
vertex. Since at the LHC the searches will rely heavily on producing the Higgs boson in 
gluon-gluon fusion, and since the mass determination will rely in part on the decays into 
gluon pairs, such a scenario may present experimental difficulties. 

• In the small-a e ff scenario the couplings governing the decays h^ bb and h^ t + t~ 
are suppressed with respect to their Standard Model values by a factor — sma e s/ cos (3 
(ojcfr is the effective mixing angle of the neutral CP-even Higgs sector including radiative 
corrections). The suppression occurs mainly for large tan/5 and moderate tua- 

2.2 CP-violating scenarios 

In CP-violating MSSM scenarios the three neutral Higgs mass eigenstates Hi (i = 1, 2, 3) do not 
have well defined CP quantum numbers. Each of them can thus be produced by Higgsstrahlung 
(e + e - ^ HiL) via the CP-even field component and in pairs (e + e - -^ HiHj (i ^ j)). The 
relative rates depend on the choice of the parameters describing the CP-even/odd mixing. 

Experimentally, the CP-violating scenarios are more challenging than the CP-conserving 
scenarios. For a wide range of model parameters, the coupling of the lightest Higgs boson 
Tii to the Z boson may be suppressed. Furthermore, the second- and third-lightest 7i 2 and 
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H.3 bosons may both have masses close to or beyond the kinematic reach of LEP. Also, in 
CP- violating scenarios, the decays to the main "discovery channels", Ti\— > bb, 7i 2 ^ bb and 
TC 2 — > T~L\H\—* bbbb 4 , may have lower branching ratios. One therefore anticipates less search 
sensitivity in the CP-violating scenarios than in the CP-conserving scenarios. An example 
illustrating this situation is given in Table El 



Parameters 


FeynHiggs 


CPH 


H+ (GeY/c 2 ) 


129.0 


129.0 


tan (3 


5.0 


5.0 


m Hl (GeV/c 2 ) 


38.1 


33.4 


m H2 (GeV/c 2 ) 


105.4 


102.4 


<r(WiZ-> bbZ) (pb) 


0.0051 


0.0019 


a(H 2 Z^ bbZ) (pb) 


0.0156 


0.0197 


a(H 2 Z^ HxHxL^ bbbbZ) (pb) 


0.0866 


0.0978 


oXuxU^ bbbb) (pb) 


0.0066 


0.0094 



Table 2: A typical parameter set which is difficult to address by the present searches. The 
results of the two calculations, FeynHiggs and CPH, are given for a centre-of-mass energy of 
206 GeV. The main input parameters are listed in the first two lines; all other input parameters 
correspond to the CPX benchmark scenario and are listed in the last column of Table The 
output masses m-^ 1 , m-^ 2 and the relevant topological cross-sections are listed below the second 
horizontal line. 



The cross-sections for Higgsstrahlung and pair production are given by [9] 

^Z = 9 Ul ZZ °WL I 5 ) 

2 \ SM (n\ 

O-Hi-Hi = 9-HiHjZ X °WL (0) 

(in the expression for A, Eq. 3, the indices h and A have to be replaced by Tii and TCj). The 
couplings 

gn.zz = cos (30u + sin/30 2 i (7) 
= 3i {cos (30 2j - sin (30^) - 3j (cos (30 2l - sin (30u) (8) 
obey the complementarity relation 

£<4zz = l (9) 
1=1 

4 Regarding the decay properties, the CP-violating scenarios maintain a certain similarity to the CP- 
conserving scenarios although the branching ratios are, in general, different. The lightest mass eigenstate 
Tlx predominantly decays to bb if allowed by kinematics, with a small fraction decaying to t + t~ and cc. The 
second-lightest Higgs boson H2 may decay to TiiHi when allowed by kinematics; otherwise it decays preferen- 
tially to bb. 
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9n k ZZ — ^ijkg-HiHjZ 



(10) 



where is the usual Levi-Civita symbol. 

In CP- violating scenarios, the orthogonal matrix Oy (i,j = 1,2,3) relating the weak CP 
eigenstates to the mass eigenstates has non-vanishing off-diagonal elements. These elements, 
giving rise to CP-even/odd mixing, are proportional to 



with v — \Jvi + v\. Substantial deviations from the CP-conserving scenarios are thus expected 
for small Msusy anc ^ l ar g e I m (/zv4), which are obtained if the CP- violating phase arg(A) takes 
values close to 90°. Furthermore, the effects from CP violation strongly depend on the precise 
value of the top quark mass [15]. 

The parameters of the benchmark model CPX have been chosen [18] to maximise the phe- 
nomenological differences with respect to the CP-conserving scenarios. Constraints from mea- 
surements of the electron and neutron electric dipole moments [26] were also taken into account. 
The basic set of parameters is listed in the last column of Table ^ Note that the scan of mn± 
started at 4 GeV/c 2 but values less than about 100 GeV/c 2 give unphysical results and are 
thus considered as theoretically inaccessible. 

The parameters which follow have been varied one-by-one while all the other parameters 
were kept at their standard CPX value. 

• Top quark mass: m t = 169.3,174.3,179.3 and 183.0 GeV/c 2 , embracing the current 
experimental value, m t = 172.7 ± 2.9 GeV/c 2 [15]. 

• The CP-violating phases: arg(A) = arg(mg) = 0°, 30°, 60°, 90° (CPX value), 135° and 
180° (the values 0° and 180° correspond to CP-conserving limits). 

• The Higgs mass parameter: /x = 0.5, 1.0, 2.0 (CPX value) and 4.0 TeV. 

• The SUSY-breaking scale: M SU sy = 0.5 TeV (CPX value) and 1.0 TeV. The proposal of 
the CPX scenario [18] predicts a weak dependence on Msusy if the relations \A\ = \mg\ = 
/i/2 = 2Msusy are preserved. This behaviour is examined by studying a model where 
-^susy is increased from 0.5 TeV to 1 TeV and the values of A, mg and \i are scaled to 
2000 GeV, 2000 GeV and 4000 GeV, respectively. 



3 Experimental searches 

The searches carried out by the four LEP collaborations are based on e + e~ collision data 
which span a large range of centre-of-mass energies, from 91 GeV to 209 GeV. The searches 
include the Higgsstrahlung and pair production processes, ensuring by their complementarity 



m\ Im(fiA) 



(11) 
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a high sensitivity over the accessible MSSM parameter space. It is important to note that 
the kinematic properties of the signal processes are to a large extent independent of the CP 
composition of the Higgs bosons. This implies that the same topological searches can be applied 
to study the CP-conserving and CP-violating scenarios. For Higgsstrahhmg this is natural since 
only the CP-even components of the Higgs fields couple to the Z boson. In pair production 
involving CP-even and CP-odd field components, the similarity of the kinematic properties 
(e.g., angular distributions) arises from the scalar nature of the Higgs bosons. Small differences 
may occur from spin-spin correlations between final-state particles but these were found to have 
no noticeable effect on the signal detection efficiencies. We therefore adopt in the following a 
common notation for the CP-conserving and CP-violating processes in which Hi (i = 1,2,3) 
designate three generic neutral Higgs bosons of increasing mass, with undefined CP properties; 
in the CP-conserving limit (arg(A) = arg(mg) = 0°), these become the CP eigenstates h, A, H 
(the correspondence depends on the mass hierarchy). 

In each of the four LEP experiments, the data analysis is done in several steps. A pres- 
election is applied to reduce some of the largest backgrounds, in particular, from two-photon 
processes. The remaining background, mainly from production of fermion pairs and WW or 
ZZ (possibly accompanied by photon or gluon radiation), is further reduced by more selective 
cuts or by applying multivariate techniques such as likelihood analyses and neural networks. 
The identification of b-quarks in the decay of the Higgs bosons plays an important role in 
the discrimination between signal and background, as does the kinematic reconstruction of the 
Higgs boson masses. The detailed implementation of these analyses, as well as the data samples 
used by the four collaborations, are described in the individual publications. A full catalog of 
the searches provided by the four LEP collaborations for this combination, with corresponding 
references to the detailed descriptions, is given in Appendix A. 

3.1 Search topologies 

Searches have been carried out for the two main signal processes, the Higgsstrahlung process 
e + e~— > Tii Z (which also apply in some cases to e + e~^ T-^Z) and the pair production process 

e+e"^ n 2 n 1 . 

(a) Considering first the Higsstrahlung process e + e~— > TiiZ, the principal signal topologies are 
those used in the search for the Standard Model Higgs boson at LEP [3], namely: 

• the four-jet topology, (7~li—> bb)(Z^ qq), in which the invariant mass of two jets is close 
to the Z boson mass M z while the other two jets contain b-flavour; 

• the missing energy topology, {Ti\— > bb, r + r~)(Z-^ uu), in which the event consists of 
two b-jets or identified tau decays and substantial missing momentum and missing mass, 
compatible with M z ; 

• the leptonic final states, (T~ti— > bb)(Z— > e + e~, /z + /i~), in which the invariant mass of the 
two leptons is close to M z ; 
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• the final states with tau-leptons, (Hi— > t + t )(Z— > qq) and (Hi— > bb, r + r )(Z— > r + r ), 
in which either the t + t~ or the qq pair has an invariant mass close to M z . 

Most of these signatures are relevant for Higgs boson masses above the bb threshold and rely 
on the identification of b-quarks in the final state. Searches for lighter Higgs bosons, listed in 
Appendix A, use signatures which are described in the specific publications. In some regions 
of the MSSM parameter space, the 7Yi— > bb decay may be suppressed while decays into other 
quark flavours or gluon pairs are favoured. The above searches are therefore complemented or 
replaced 5 by flavour-independent searches for (Hi— > qq)Z in which there is no requirement on 
the quark-flavour of the jets. Finally, the searches for Higgsstrahlung also include the Higgs 
cascade decay e + e~-^ H{L^> (Tli Hi)Z, giving rise to a new class of event topologies. These 
processes may play an important role in those regions of the parameter space where they are 
allowed by kinematics. 

(b) In the case of the pair production process, e + e~-^ H2H1, the principal signal topologies 
at LEP are: 

• the four-b final state (H 2 — ► bb)(Hi— > bb); 

• the mixed final states (H.2~^ t + t~) (Hi—> bb) and (H 2 ^> bb)(7ii— > t + t~); 

• the four-tau final state (Ti.2—> t + t~)(H\^ t + t~). 

The Higgs cascade decay, e + e~-^ H 2 H\^ (HiHi)Hi, gives rise to event topologies ranging 
from six b-jets to six tau-leptons. Most of these searches are relevant for Higgs boson masses 
above the r + r~ threshold. Similarly to the Higgsstrahlung case, the above searches for pair 
production are complemented or replaced, whenever more efficient, by flavour-independent 
searches. 



3.2 Additional experimental constraints 

If the combination of the above searches is not sufficiently sensitive for excluding a given model 
point, additional constraints are applied; these are listed below. 

• Constraint from the measured decay width of the Z boson, T z , and its possible deviation, 
Ar z , from the Standard Model prediction. The model point is regarded as excluded if 
the following relation between the relevant cross-sections is found to be true: 

Ar 

E^N + E^W >-^-*r(mz), (12) 

1 7 

5 The replacement is necessary whenever the overlap in terms of selected events is important, in order to 
avoid double-counting. 
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where ATz = 2.0 MeV [27] stands for the 95% CL upper bound on the possible additional 
decay width of the Z boson, beyond the Standard Model prediction, and a z ot is the Z pole 
cross-section. 

• Constraint from a decay mode independent search for e + e~— > 7i\L [28]. The model point 
is regarded as excluded if the condition 

<7„ iZ > k(m ni ) ■ a™ (13) 

is fulfilled, where k(m Hi ) is a mass-dependent factor which scales the Standard Model 
Higgs production cross-section to the value that is excluded at the 95% CL. 

• Constraint from a search for light Higgs bosons produced by the Yukawa process 6 . The 
model point is regarded as excluded if the predicted Yukawa enhancement factor £(m Wl ), 
defined in [29], is excluded by this search. To be conservative, the weaker of the two 
enhancement factors, for CP-even and CP-odd couplings, is used. 

These additional constraints are particularly useful at small m-^ 1 and m^ 2 , below the bb thresh- 
old. 



3.3 Statistical combination of search channels 

The statistical method by which the topological searches are combined is described in Refs. [3, 
30]. 

After selection, the combined data configuration (distribution of all selected events in several 
discriminating variables) is compared in a frequentist approach to a large number of simulated 
configurations generated separately for two hypotheses: the background (b) hypothesis and the 
signal-plus-background (s + b) hypothesis. The ratio 

Q = C s+b /C b (14) 

of the corresponding likelihoods is used as the test statistic. The predicted, normalised, 
distributions of Q (probability density functions) are integrated to obtain the p- values [31] 
1 - CL b = 1 - V b (Q < ^observed) and CL s+b = V s+b (Q < ^observed); these measure the com- 
patibility of the observed data configuration with the two hypotheses. Here V b and V s+b are 
the probabilities for a single experiment to obtain a value of Q smaller than or equal to the 
observed value, given the background or the signal-plus-background hypothesis. More details 
can be found in Ref. [3]. 

Systematic errors are incorporated in the calculation of the likelihoods by randomly varying 
the signal and background estimates in each channel 7 according to Gaussian error distributions 

6 Note that, in the case of DELPHI, the Yukawa channels are not used as external constraints but are 
combined with the other search channels. 

7 The word "channel" designates any subset of the data in which a search has been carried out. These subsets 
may correspond to specific final-state topologies, to data sets collected at different centre-of-mass energies or to 
the subsets of data collected by different experiments. 
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and widths corresponding to the systematic errors. For a given source of uncertainty, correla- 
tions are addressed by applying these random variations simultaneously to all those channels 
for which the source of uncertainty is relevant. Errors which are correlated among the experi- 
ments arise mainly from using the same Monte Carlo generators and cross-section calculations 
for the signal and background processes. The uncorrelated errors arise mainly from the limited 
statistics of the simulated background event samples. 

In a purely frequentist approach, the exclusion limit is computed from the confidence CL s+b 
for the signal-plus-background hypothesis: a signal is regarded as excluded at the 95% CL, for 
example, if an observation is made such that CL s+b is lower than 0.05. However, this procedure 
may lead to the undesired situation in which a large downward fluctuation of the background 
would exclude a signal hypothesis for which the experiment has no sensitivity since the expected 
signal rate is too small. This problem is avoided by using the ratio 

CL S = CL s+b /CL b (15) 

instead of CL s+b . We adopt this quantity for setting exclusion limits and consider a given model 
to be excluded at the 95% CL if the corresponding value of CL S is less than 0.05. Since CL b 
is a positive number less than one, CL S is always larger than CL s+b and the limits obtained in 
this way are therefore conservative. 

3.4 Comparisons of the data with the expected background 

The distribution of the p- value 1 — CL b over the parameter space covered by the searches 
provides a convenient way of studying the agreement between the data and the expected back- 
ground and of discussing the statistical significance of any local excess in the data. While a 
purely background-like behaviour 8 would yield p-values close to 0.5, much smaller values are 
expected in the case of a signal-like excess. For example, a local excess of three or five standard 
deviations would give rise to a p-value 1 — CL b of 2.7 x 10~ 3 or 5.7 x 10~ 7 , respectively. 

One has to be careful, however, when interpreting these numbers as probabilities for lo- 
cal excesses occurring over the extended domains covered by the searches. For example, the 
probability for a fluctuation of three standard deviations to occur anywhere in the parameter 
space is much larger than the number 2.7 x 10~ 3 just quoted. A multiplication factor has 
to be applied to the probability 1 — CL b which reflects the number of independent "bins" of 
the parameter space; this factor can be estimated from the total size of the parameter space 
and the experimental resolutions. For example, the searches for the Higgsstrahlung process 
e + e~— > HiZ, covering the range 0<m^ 1 <120 GeV/c 2 with a mass resolution Am^ x of about 
3 GeV/c 2 , would yield about twenty fairly independent mass-bins of width 2Am-^ 1 ; hence, a 
multiplication factor of about twenty. Much bigger multiplication factors are expected in the 
searches for the pair production process e + e~— > 7i 2 7^i with two independent search parameters 
(masses) . 

8 Single, background-like, experiments have values of 1 — CL b uniformly distributed between zero and one. 
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These simple considerations do not take into account, for example, possible correlations 
from resolution tails extending over several adjacent bins or correlations between different 
searches sharing candidate events. A more elaborate evaluation of the multiplication factor has 
therefore been performed. A large number of background experiments was simulated, covering 
the whole parameter space, using realistic resolution functions and taking correlations into 
account. From these random experiments, the probability to obtain 1 — CL^ smaller than a 
given value, anywhere in the parameter space of a given scenario, has been determined (the 
m^-max scenario was taken for this study). A scale factor of at least 60 was obtained in this 
manner. According to this estimate, the probability of observing a background fluctuation of 
three standard deviations anywhere in the parameter space of a given scenario (e.g., m^-max) 
can be 16% or more. Also, to observe two fluctuations with two standard deviations turns out 
to be more likely than to observe only one. 

Figure^shows the distribution of the p- value 1— CLj,, determined from the present combined 
searches, for the CP-conserving benchmark scenario m^-max and the CP- violating scenario 
CPX. Over the largest part of the parameter space, the local excesses are smaller than two 
standard deviations. In the m^-max scenario, the lowest value, 1 — CLb = 1.3 x 10~ 2 , lies 
within the vertical band at around 100 GeV/c 2 and corresponds to 2.5 standard deviations. 
This excess, and a less significant excess at about 115 GeV/c 2 , come from the Higgsstrahlung 
search; both are discussed in Ref. [3] in the context of the search for the Standard Model 
Higgs boson. In the CPX scenario, one observes two small regions at rs 35-40 GeV/c 2 , 
m-}i 2 ~ 105 GeV/c 2 and tan/3 ~ 10, where the significance exceeds three standard deviations; 
they arise from the search for the pair production process. 

The exact position and size of these fluctuations may vary from one scenario to the other. 
In Tables H3 and 0] we list the parameters of the most significant excesses for all CP-conserving 
and CP-violating benchmark scenarios considered in this paper. The largest fluctuation of all 
has a significance of 3.5 standard deviations; its probability is estimated as 3.6% at least, when 
the scale factor of 60 or more is applied. 

From these studies one can conclude that there is a reasonable agreement between the data 
and the simulated background, with no compelling evidence for a Higgs boson signal, and that 
the excesses observed are compatible with random fluctuations of the background. 



4 Limits on topological cross-sections 

In this section we present upper bounds on the cross-sections for the most important final-state 
topologies expected from the Higgsstrahlung process e + e - -^ 7i\L and the pair production 
process e + e - -^ Ji-iHx. These can be used to test a wide range of specific models. 

We define the scaling factor 

S95 — ^max/^ref > (16) 

where <r max is the largest cross-section compatible with the data, at the 95% CL, and cr rc f 
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Table 3: The most significant excesses with respect to the predicted background, for each of the 
CP-conserving benchmark scenarios. Columns 2 to 6 show the mass parameters (in GeV/c 2 ) 
and tan/3 at which the excess occurs. Column 7 gives the corresponding p-values 1 — CLf,. In 
the last column, the signihcances of the excesses, in standard deviations, are listed. 



is a reference cross-section. For the topologies motivated by Higgsstrahlung, <7 re f is taken to 
be the Standard Model Higgs production cross-section; for final states motivated by the pair 
production process, cr re f is taken to be the MSSM Higgs production cross-section of Eq. 2 with 
the MSSM suppression factor set to 1. Numerical values for the cross-section limits are listed 
in Appendix B. 

Figure El shows the upper bound S95 for final states motivated by the Higgsstrahlung process 
e + e~— > TiiL (the figure is reproduced from Ref. [3]). In part (a), the Higgs boson is assumed to 
decay into fermions and bosons with branching ratios as given by the Standard Model. Contri- 
butions from the fusion processes WW^ Hi and ZZ^ TCi, according to the Standard Model, 
corrected for initial-state radiation, are assumed to scale with energy like the Higgsstrahlung 
process. In part (b) it is assumed that the Higgs boson decays exclusively to bb and in part (c) 
exclusively to t + t~. Besides representing bounds on topological cross-sections, this figure also 
illustrates the overall agreement between the data and the expected background from Standard 
Model processes. The largest deviations observed barely exceed two standard deviations. 

Figure El shows contours of S95 for the cascade process e + e"-^ 7Y2Z— > (TC\TCi)Z, projected 
onto the (m^ 2 , m^J plane, assuming that the TC2 boson decays exclusively to TC1TC1. In part 
(a) it is assumed that the Hi boson decays exclusively to bb and in part (b) exclusively to 
t + t~ . In part (c), as an example, an equal mixture of Hi— > bb and Hi— > t + t~ is assumed, 
which implies 25% bbbbZ, 25% t + t~t + t~Z and 50% bbr + T~Z final states. The sensitivity of 
the bbbbZ channel starts at the bb threshold and extends almost to the kinematic limit. In the 
t + t~t + t~Z channel the sensitivity is altogether weaker (the discontinuities reveal the limited 
and inhomogeneous mass coverage of the four experiments in this channel). 
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Table 4: The most significant excesses with respect to the predicted background in the CP- 
violating benchmark scenario CPX and its variants. The first column indicates either the 
CPX scenario or the parameter value which differs from the standard CPX set listed in the 
last column of Table Columns 2 to 6 show the mass parameters (in GeV/c 2 ) and tan/? 
at which the excesses occur (the more conservative of the CPH and FeynHiggs calculations is 
used). Columns 7 and 8 give the corresponding p-values, 1 — CL b , using in turn the CPH and 
FeynHiggs codes (note the overall agreement of the two calculations in this respect). In the 
last column, the significances of the excesses, in standard deviations, are listed. 



Figure |U shows S95 for final states motivated by the pair-production process e + e~— > H2H1, 
for the particular case where the masses m-^ 2 and m-H\ are approximately equal. Such is the 
case, for example, in the CP-conserving MSSM scenario m^-max for tan/3 larger than about 10 
and small m^ 2 (= tua). In part (a), the 7^2 and H\ decay branching ratios correspond to the 
m^-max benchmark scenario with tan/3 = 10 (see the caption for the exact values); in part (b), 
both 7i 2 and Hi are assumed to decay exclusively to bb; in part (c), one Higgs boson is assumed 
to decay exclusively to bb while the other exclusively to t + t"; in part (d), H2 and Hi are both 
assumed to decay exclusively to t + t - . At low masses, the exclusion limits are completed using 
the constraint from the measured decay width of the Z boson (see Section 3.2). This figure also 
illustrates the overall agreement between the data and the expected background from Standard 
Model processes since the largest deviations are within two standard deviations. 

Figure El shows contours of S95 for final states motivated by the process e + e~^ H2H1, 
projected onto the {mn 2 , mm) plane. In part (a), both Higgs bosons are assumed to decay 
exclusively to bb and in part (b) exclusively to t + t~. In parts (c) / (d), the H2 / Hi boson 
is assumed to decay exclusively to bb while the other boson is assumed to decay exclusively to 



15 



Figure |H1 shows contours of S 95 for the cascade process e + e~— > Ti^Hi^- (Hi'Hi)'Hi, projected 
onto the (rnn 2 , m-Hi) plane, assuming that the H2 boson decays exclusively to H.{K\. In part 
(a), the Hi boson is assumed to decay exclusively to bb and in part (b) exclusively to t + t~. 
In part (c), as an example, an equal mixture of Hi— > bb and Ti\— > t + t~ is assumed, which 
implies 12.5% bbbbbb, 37.5% bbbbY+r~, 37.5% bbr+r-r+r- and 12.5% t + t-t + t-t + t- final 
states. 

A word of caution is in place concerning the correlations which exist between some of the 
above cross-section limits which arise from overlapping candidates in the corresponding selec- 
tions. Such correlations are present, for example, between b-tagged and flavour-independent 
searches of a given experiment or between searches addressing direct decays (e.g., TL{L-^ bbbb) 
and cascade decays (e.g., (7^2^ TCiTCi)Z—> bbbbbb); they may be a source of problems if several 
of the cross-section limits are used in conjunction to test a given model. Note, however, that 
these correlations are properly taken into account in the model interpretations which follow. 



5 Results interpreted in CP-conserving MSSM scenarios 

In this section, the search results are interpreted in the CP-conserving benchmark scenarios 
presented in Section 2.1. The exclusion limits, which are shown in the figures below at the 95% 
CL and the 99.7% CL, are obtained from the values of CL S (see Eq. 15), for an assumed top 
quark mass of m t = 174.3 GeV/c 2 . The exclusion limits are presented in four projections of 
the MSSM parameter space. The limits expected on the basis of Monte Carlo simulations with 
no signal, at the 95% CL, are also indicated. The exact mass bounds and exclusions for tan f3 
are listed in Table El for four values of m t . 

The exclusions for the m^-max benchmark scenario are shown in Figure [7| In the region 
with tan (3 less than about five, the exclusion is provided mainly by the Higgsstrahlung process, 
giving a lower bound of about 114 GeV/c 2 for m h . At high tan/3, the pair production process 
is most useful, providing limits in the vicinity of 93 GeV/c 2 for both and toa- For in the 
vicinity of 100 GeV/c 2 , one observes a deviation between the expected and the experimental 
exclusions. This deviation, which is also present in other CP-conserving scenarios, is due to the 
excess in the Higgsstrahlung channel which was discussed in Ref. [3] and gives rise to the vertical 
bands in Figures ^ (a) and (b). Note that the mass bounds obtained are largely insensitive to 
the top quark mass. 

The data also exclude certain domains of tan/3. This is best illustrated in the (mh, tan/3) 
projection (plot (b)) where the upper boundary of the parameter space along is indicated 
for four values of m t ; the intersections of these boundaries with the experimental exclusions 
define the regions of tan /3 which are excluded. The exclusion in tan (3, as a function of the 
assumed top quark mass, is summarised in Figure |H1 for m t larger than about 181.5 GeV/c 2 , 
no 95% CL limit on tan (3 can be set in this scenario. 
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One should be aware that the upper boundary of the parameter space along also depends 
moderately on the choice of Msusy- For example, changing Msusy from 1 TeV to 2 TeV would 
broaden the parameter space by about 2 GeV/c 2 along mh, with corresponding effects on the 
exclusions in tan /?. This observation holds for all CP-conserving scenarios which follow. 

Figures 19*1 and ITUl show the same set of plots for the two variants, (a) and (b), of the m^-max 
scenario introduced in Section 2.1.1. The change of the sign of the Higgs mass parameter fi or 
of the mixing parameter X t barely affect the mass limits; however, sizable differences occur in 
the exclusions of tan/? (see Table EJ). For example, in variant (b), a small domain of tan/? is 
excluded even for m t = 183.0 GeV/c 2 , which is not the case in the standard m^-max scenario 
and its variant (a). Note, in Figure |""J the small domains at between 60 and 75 GeV/c 2 , 
small toa and tan/?<0.9 which are excluded at the 95% CL but not at the 99.7% CL. 

The exclusions for the CP-conserving no-mixing benchmark scenario are shown in Figure ITT1 
In this scenario, the theoretical boundaries of the parameter space are more restricted than in 
the m^-max scenario. As a consequence, large domains of tan j3 are excluded for all the top 
quark masses considered. Note the relatively strong variation of the exclusion limits with mt in 
this scenario (see Table EJ), which is caused by the proximity of the experimental lower bound 
of mh from the Higgsstrahlung searches and the theoretical upper bound of m^. 

An interesting feature of this scenario is that, for larger than about 100 GeV/c 2 and 
large tan /?, the heavy scalar boson H is within kinematic reach. Moreover, the cross-section 
for the process e + e~^ HZ is increasing with tan/3, resulting in an improved search sensitivity; 
this explains the nearly circular shape of the expected limit in Figure ITT1 (b). 

Note the small domain at between 75 and 80 GeV/c 2 , small m\ and tan/? < 0.7, barely 
perceptible in the plots, which is not excluded in this scenario at 95% CL (this domain is 
excluded for m t = 169.3 GeV/c 2 ). The branching ratio for h— > bb is small and the decay 
h— > AA is dominant in this region. The A boson, with mass below the t + t~ threshold, may 
decay to final states which are not sufficiently covered by the present searches. For this reason, 
the mass limits given in TableEJfor this scenario and for mt larger than 169.3 GeV/c 2 are valid 
only for tan/? > 0.7. Conversely, for m t larger than 169.3 GeV/c 2 , the quoted exclusion of 
tan/? is valid only for rriA larger than about 3 GeV/c 2 . 

Figure shows the exclusion plots for the (a) variant of the no-mixing scenario introduced 
in Section 2.1.2. The change of sign of the Higgs mass parameter [i and the increase of the weak 
SUSY-breaking scale from 1 TeV to 2 TeV affect only the theoretical bounds of the parameter 
space but barely change the mass limits, except for m t =169.3 GeV/c 2 . There are moderate 
changes though in the exclusions of tan /?. In the hatched domain (tan /?<0.7), the contributions 
from top and stop quark loops to the radiative corrections are large and uncertain; hence, no 
exclusions can be claimed there. 

The exclusions for the large-n benchmark scenario are shown in Figure [j"""J As mentioned 
in Section 2.3, this scenario was constructed to test the sensitivity of LEP to MSSM scenarios 
which may be a priori difficult to handle experimentally since the Higgs boson decays to bb 
are largely suppressed. It turns out that the flavour-independent and decay-mode-independent 
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searches are sufficiently powerful to exclude all such situations at 95% CL, for top quark masses 
up to 174.3 GeV/c 2 . There remains a thin strip at tan/? larger than about 10 and running from 
niA of about 100 to about 200 GeV/c 2 , which is excluded at the 95% CL but not at 99.7% CL 
because the suppression of the bb channel is particularly strong in that region. This strip is 
found to grow with increasing m t and becomes gradually non-excluded at the 95% CL. Other 
small, weakly excluded, regions are located at m h « 60 GeV/c 2 and small m^, and along the 
mh « 771a "diagonal" of plot (a). 

Similar plots are shown in Figures mi and [TBI for the gluophobic and small-a e ff scenarios 
defined in Section 2.1.3. These scenarios were designed to test situations which can be prob- 
lematic at the Tevatron and LHC colliders. In both cases, large domains of the parameter space 
are excluded by the LEP searches. 



6 Results interpreted in CP-violating MSSM scenarios 

In this section, the search results are interpreted in the CP-violating benchmark scenario CPX 
presented in Section 2.2, and in some variants of CPX where the basic model parameters are 
varied one-by-one. Note that in these scenarios mu 3 is always larger than 120 GeV/c 2 , except 
where the CP-violating phases arg(A) = arg(mg) are put to 0° or 180°. 

The experimental exclusions for the CPX benchmark scenario are shown in Figure EE in 
four projections. For large m^ 2 , the TCi is almost completely CP-even; in this case the limit 
on m?^ is close to 114 GeV/c 2 , the limit obtained for the Standard Model Higgs boson [3]. 
For example, for m-H 2 larger than 133 GeV/c 2 , one can quote a lower bound of 113 GeV/c 2 for 
m-Hx- Large CP-odd admixtures to Tii occur, however, for smaller mn 2 , giving rise to domains 
at lower which are not excluded. 

The exclusion is particularly weak for tan/3 between about 3.5 and 10. Here, the signal 
is spread over several channels arising from the Higgsstrahlung and pair-production processes, 
including the 7iiH\ cascade decays, which give rise to complex final states with six 

jets. The parameter set of Table El is a typical example of this situation. This is illustrated 
in Figure El where the main final-state cross-sections are plotted as a function of tan (5 (the 
FeynHiggs calculation is used). In general, these signal contributions cannot be added up 
statistically because of a large overlap in the selected events; hence, a relatively low overall 
detection efficiency is expected. Moreover, one of the experiments presents a local excess of 
about two standard deviations in this domain of tan/3 and for of about 45 GeV/c 2 [10], 
which lowers the exclusion power below the expectation. Nonetheless, the region defined by 
m?^ < 114 GeV/c 2 and tan/3 < 3.0 is excluded by the data (see Figure ITBI (b)) and a 95% CL 
lower bound of 2.9 can be set on tan/3 in this scenario. 

Figure HH1 illustrates the exclusions in the (m^, tan/3) projection, using the CPH calculation 
(part (a)) and the FeynHiggs calculation (part (b)). Differences occur mainly at large tan/? 
where the FeynHiggs calculation predicts a larger Higgsstrahlung cross-section and hence a bet- 
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ter search sensitivity than the CPH calculation. In parts (a) and (b) of the figure, one observes 
two distinct domains at moderate tan (3, with m Wl <15 GeV/c 2 and 30 GeV/c 2 <m Wl <55 GeV/c 2 , 
which are not excluded at the 95% CL. The values of 1 — CLb indicate that these domains are 
excluded, respectively, at the 55% CL and 77% CL using the CPH calculation, and at the 50% 
CL and 66% CL, respectively, using the FeynHiggs calculation. A third domain appears in 
part (b) at higher (where the CPH calculation indicates no exclusion power at all); this 
domain is excluded at the 42% CL using FeynHiggs. 

As explained in Section 2, neither of the two approaches, CPH or FeynHiggs, are preferred 
on theoretical grounds. For this reason, part (c) of this figure was obtained by choosing in each 
scan point of the parameter space the more conservative of the two approaches, i.e., the one for 
which the less significant exclusion is observed. The same procedure was adopted in Figure 
and in all the figures which follow. 

The significant impact of the top mass on the CP-violating effects, indicated by Eq. 11, 
is illustrated in Figure IT^l where the (m^, tan/?) projection is shown for four values of 
m t . With increasing m t , one observes a reduction of the exclusion power, especially in the 
region of tan/5 between 3.5 and 10. No lower bound on can be quoted in this do- 

main. In plot (a) (for m t = 169.3 GeV/c 2 ), the two domains with < 15 GeV/c 2 and 
30 GeV/c 2 < m Hl < 55 GeV/c 2 are excluded at the 60% CL and 88% CL, respectively. 

Figure I2*U1 illustrates the exclusion in the (m^, tan/?) plane as a function of the CP- 
violating phases, axg(A) = arg(m g ), which are varied together. For phase angles close to 
0°, the experimental exclusions are similar to those in the CP-conserving scenarios (see, for 
example, Figure [3 but note the differences in the allowed parameter space). Sizable differences 
are observed for larger phase angles, especially for arg(A) = arg(mg) = 90° (the CPX value). 
At arg(A) = arg(mg) = 180° (another CP-conserving scenario), the allowed parameter space is 
excluded almost completely. Note however that in the hatched region, with tan (3 greater than 
about 12, the calculation of the bottom- Yukawa coupling has large theoretical uncertainties; 
hence no exclusion can be claimed in this domain. 

In Figure the value of the Higgs mass parameter /i is varied from 500 GeV through 
1000 GeV and 2000 GeV (the CPX value) to 4000 GeV. At small values, the CP-violating 
effects are small (see Eq. 11) and the exclusion power is strong (as in the CP-conserving case). 
For fi larger than 2000 GeV and large tan j3, the FeynHiggs and CPH calculations both provide 
bottom- Yukawa coupling in the non-perturbative regime, giving rise to negative values for the 
square of and to other unphysical results. For /i < 2000 GeV this regime sets in only at 
tan j3 larger than 40 whereas for /i = 4000 GeV this situation already occurs at tan f3 abowe 
20. Hence, in Figure I2~T1 (d), the hatched domain should not be considered as being integrally 
part of the allowed parameter space. 

Figure E21 illustrates the dependence on the soft SUSY-breaking scale parameter, M SU sy, 
which is increased from the CPX value of 500 GeV in part (a) to 1000 GeV in part (b). This 
decreases the CP- violating effects (see Eq. 11) and leads to a larger exclusion. The "scaling" 
behaviour mentioned in Section 2.3, namely the relative insensitivity of the exclusions to changes 
in M S usy as long as the relations |A,t>| = \mg\ = A*/2 = 2M SU sy are preserved, is qualitatively 
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confirmed by comparing parts (a) and (c) of the figure. 



7 Summary 

The searches for neutral Higgs bosons described in this paper are based on the data collected 
by the four LEP collaborations, ALEPH, DELPHI, L3 and OPAL, which were statistically 
combined by the LEP Working Group for Higgs Boson Searches. The data samples include 
those collected during the LEP 2 phase at e + e~ centre-of-mass energies up to 209 GeV; two 
experiments also provided LEP 1 data, at energies in the vicinity of the Z boson resonance. 
The searches address a large number of final-state topologies arising from the Higgsstrahlung 
process e + e~^ 7i\L and from the pair production process e + e - -^ IH.2H1. The combined LEP 
data do not reveal any excess of events which would indicate the production of Higgs bosons. 
The differences with respect to the background predictions are compatible with statistical 
fluctuations of the background. 

From these results, upper bounds are derived for the cross-sections of a number of Higgs- 
like event topologies. These upper bounds cover a wide range of Higgs boson masses and are 
typically well below the cross- sect ions predicted within the MSSM framework; these limits can 
be used to constrain a large number of theoretical models. 

The combined search results are used to test several MSSM scenarios which include CP- 
conserving and CP- violating benchmark models. These models are motivated mainly by physics 
arguments but some of them are constructed to test specific situations where the detection of 
Higgs bosons at the Tevatron and LHC colliders might present experimental difficulties. It is 
found that in all these scenarios the searches conducted at LEP exclude sizable domains of the 
theoretically allowed parameter space. 

In the CP-conserving case, lower bounds can be set on the masses of neutral Higgs bosons 
and the value of tan/3 can be restricted. Taking, for example, the CP-conserving scenario 
m^-max and a top quark mass of 174.3 GeV/c 2 , values of and less than 92.8 GeV/c 2 
and 93.4 GeV/c 2 , respectively, are excluded at the 95% CL. In the same scenario, values of 
tan/? between 0.7 and 2.0 are excluded, but this range depends considerably on the assumed 
top quark mass and may also depend on M S usy- 

In the CP-violating benchmark scenario CPX and the variants which have been studied, the 
combined LEP data show large domains which are not excluded, down to the lowest mass values; 
hence, no absolute limits can be set for the Higgs boson masses. The excluded domains vary 
considerably with the precise value of the top quark mass and the MSSM model parameters. 
For example, in the CPX scenario with standard parameters and m t = 174.3 GeV/c 2 , tan/3 
can be restricted to values larger than 2.9 at the 95% CL. 
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Benchmark m t (GeV/c 2 ) 
scenario 


m h (GeV/c 2 ) 


m A (GeV/c 2 ) 


Exclusions of 
tan (3 


m^-max 169.3 
174.3 
179.3 
183.0 


92.9 (94.8) 

92.8 (94.9) 

92.9 (94.8) 
92.8 (94.8) 


93.4 (95.1) 
93.4 (95.2) 

93.4 (95.2) 

93.5 (95.2) 


0.6-2.6 (0.6-2.7) 
0.7-2.0 (0.7-2.1) 
0.9-1.5 (0.9-1.6) 
no excl. (no excl.) 


uih-max 169.3 
(a) 174.3 
179.3 
183.0 


92.7 (94.9) 
92.7 (94.8) 

92.6 (94.8) 

92.7 (94.8) 


93.1 (95.1) 
93.1 (95.1) 
93.1 (95.1) 
93.1 (95.1) 


0.7-2.1 (0.7-2.2) 
0.7-2.1 (0.7-2.2) 
0.9-1.6 (0.8-1.7) 
no excl. (no excl.) 


mh-max 169.3 
(b) 174.3 
179.3 
183.0 


92.8 (94.8) 
92.6 (94.9) 

92.6 (94.8) 

92.7 (94.7) 


93.2 (95.2) 
93.4 (95.1) 
93.4 (95.1) 
93.4 (95.1) 


0.5-3.3 (0.5-3.5) 
0.6-2.5 (0.6-2.7) 
0.7-2.0 (0.7-2.1) 
0.8-1.7 (0.8-1.8) 


no-mixing 169.3 
174.3 
179.3 
183.0 


excl. (excl.) 
93.6 (96.0) 
93.3 (95.0) 
92.9 (95.0) 


excl. (excl.) 
93.6 (96.4) 
93.4 (95.0) 
93.1 (95.0) 


excl. (excl.) 
0.4-10.2 (0.4-19.4) 
0.4-5.5 (0.4-6.5) 
0.4-4.4 (0.4-4.9) 


no-mixing 169.3 
(a) 174.3 
179.3 
183.0 


93.2 (95.2) 
92.8 (94.9) 

92.8 (94.9) 

92.9 (94.8) 


93.4 (95.4) 
93.1 (95.1) 
93.1 (95.0) 
93.1 (95.0) 


0.7-7.1 (0.7-9.3) 
0.7-4.6 (0.7-5.1) 
0.7-3.5 (0.7-3.8) 
0.7-3.0 (0.8-3.2) 


large- u 169.3 
174.3 
179.3 
183.0 


excl. (excl.) 
excl. (excl.) 
109.2 (109.2) 
95.6 (95.6) 


excl. (excl.) 
excl. (excl.) 
225.0 (225.0) 
98.9 (98.9) 


excl. (excl.) 
excl. (excl.) 
0.7-43 (0.7-43) 
0.7-11.5 (0.7-11.5) 


gluophobic 169.3 
174.3 
179.3 
183.0 


90.6 (93.2) 
90.5 (92.3) 
90.0 (91.8) 
89.8 (91.5) 


95.7 (98.2) 
96.3 (98.0) 
96.5 (98.2) 

96.8 (98.7) 


0.4-10.3 (0.4-21.5) 
0.4-5.4 (0.4-6.4) 
0.4-3.9 (0.4-4.2) 
0.5-3.3 (0.5-3.6) 


small-a e ff 169.3 
174.3 
179.3 
183.0 


88.2 (90.0) 

87.3 (89.0) 
86.6 (88.0) 
85.6 (87.5) 


98.2 (99.6) 
98.8 (100.0) 
99.8 (100.7) 
101.0 (101.3) 


0.4-6.1 (0.4-7.4) 
0.4-4.2 (0.4-4.5) 
0.5-3.2 (0.5-3.4) 
0.6-2.7 (0.5-2.9) 



Table 5: Lower mass bounds and exclusions in tan (5, at 95% CL, obtained in the case of the 
CP-conserving MSSM benchmark scenarios, for various values of the top quark mass. In each 
case, the observed limit is followed, between parentheses, by the value expected on the basis 
of Monte Carlo simulations with no signal. In the m^-max scenario and its variant (a), there 
is no exclusion in tan/3 for m t = 183.0 GeV/c 2 or larger. The no-mixing scenario is entirely 
excluded for m t = 169.3 GeV/c 2 or smaller. In the no-mixing scenario and for m t larger than 
169.3 GeV/c 2 , the quoted mass limits are only valid for tan (3 > 0.7 and the exclusion in tan/3 
is only valid for larger than about 3 GeV/c 2 . The large- fx scenario is entirely excluded for 
m t = 174.3 GeV/c 2 or smaller. 
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Appendix A: Catalog of searches 



The searches of the four LEP collaborations which contribute to this combined analysis are 
listed in Tables El to The list is structured into two tables per experiment, one for the 
Higgsstrahlung process e + e~— >• TiiL and one for the pair production process e + e - -^ Ti^Hx- In 
each of these tables, the upper part contains the final states of the direct process and the lower 
part contains, where it applies, those of the cascade process Ti.{Hi. 

The final-state topologies are listed in the first column. In the notation adopted, Ji\ repre- 
sents the lightest and Ti.2 the second-lightest neutral Higgs boson. In the CP-conserving case, 
Tii is identified with the CP-even eigenstate h. The 7^2 is identified in most cases with the 
CP-odd eigenstate A (the cascade process TC2— > T~tiHi is identified with h— > AA). 

The symbol q indicates an arbitrary quark flavour, u, d, s, c or b. "Hadrons" include 
quarks and gluons. In the missing energy channel, in addition to the 7^iZ— > Hxvv process, the 
W fusion process l-i\v e v e (including interference) is also considered; similarly, in the leptonic 
channel, in addition to the TCiZ—> Hi£ + £~ process, the Z fusion process Hie + e~ (including 
interference) is also considered. 

The contributions based on LEP1 data (from two experiments only) can be identified by 
their value "91" in the second column which indicates the e + e~ collision energy, y/s (GeV); the 
LEP1 data used in this combination represent an integrated luminosity C of about 125 pb _1 . 
The LEP2 data span an energy range between 133 GeV and 209 GeV; they represent an inte- 
grated luminosity of about 2400 pb _1 . The integrated luminosities for the individual searches 
are listed in the third column. 

Responding to the increasing data samples and e + e~ energies, the searches were gradually 
upgraded or replaced so as to become more efficient in detecting Higgs bosons of higher masses. 
The mass ranges where the searches are relevant are listed in the next column(s). In the last 
column, references are given to the publications where the details of the searches can be found. 
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^ (GeV) 


C (pb- 1 ) 


Mass range (GeV/c z ) 


Ref. 


WiZ- (...) (...) 




mm 




(bb)(qq), (bb,cc, rr 




Oiyv) 


189 


176.2 


75 - 110 


[32] 


(any)(e + e _ , fi + n~) 






189 


176.2 


75 - 110 


[32] 


(bb)(r + r-),(r+r-)( q q) 


189 


176.2 


65 - 110 


[32] 


(bb)(qq, vv) 






192 - 202 


236.7 


60 - 120 


[33] 


(bb, r+r - , cc,gg)(e" t 


"e 


) 


192 - 202 


236.7 


60 - 120 


[33] 


(bb, r+r ,cc,gg) (r+r" 


-),(r+r-)(qq) 


192 - 202 


236.7 


60 - 120 


[33] 


(bb)(qq) 






199 - 209 


217.2 


75 - 120 


[11,34] 


(bb, r + r , cc, gg, WW) (r + r , vv) 


199 - 209 


217.2 


75 - 120 


[11,34] 


(bb,r+r~,cc,gg)(e 4 


"e 


) 


199 - 209 


217.2 


70 - 120 


[11,34] 


(bb,cc,ss,gg)(qq) 






189 


176.2 


40 - 100 


[35] 


(bb,cc,ss,gg)(z^) 






189 


176.2 


60 - 100 


[35] 


(bb,cc,ss,gg)(e+e _ , 




>-) 


189 


176.2 


60 - 115 


[32,35] 


(r+r _ )(qq) 






189 


176.2 


65 - 110 


[32] 


(bb,cc,ss,gg)(qq) 






192 - 202 


236.7 


40 - 110 


[35] 


(bb,cc,ss,gg)(z^) 






192 - 202 


236.7 


60 - 116 


[35] 


(bb,cc,ss,gg)(e+e _ , 




>") 


192 - 202 


236.7 


60 - 115 


[33, 35] 


(r+r - )(qq) 






192 - 202 


236.7 


60 - 120 


[33] 


(bb,cc,ss,gg)(qq) 






199 - 209 


217.2 


40 - 115 


[35] 


(bb,cc,ss,gg)(z/p) 






199 - 209 


217.2 


75 - 120 


[35] 


(bb,cc,ss,gg)(e+e _ , 




>-) 


199 - 209 


217.2 


70 - 120 


[11,34,35] 


(r+r" )(qq) 






199 - 209 


217.2 


60 - 120 


[11,34] 



Table 6: Summary of the ALEPH searches for the Higgsstrahlung process e+e - ^ TiiZ. The 
top part of the table lists the searches originally developed for the Standard Model Higgs boson. 
The bottom part lists flavour-independent searches where the decays of the Higgs boson into a 
quark pair of any flavour, a gluon pair or a tau pair were considered; the signal efficiencies were 
evaluated for all indicated hadronic decays of the Higgs boson. In the cases of the (T+r~) (qq) 
and leptonic channels listed in the flavour-independent part, the event selections of the Standard 
Model Higgs boson searches were used. 



29 





^(GeV) £ (pb- 1 ) 


Mass range (GeV/c 2 ) 


Ref. 


n 2 n 1 ^ (...)(...) 




(m H2 +m Hl )/2 




(bb)(bb), (r+r-)(bb), (bb)(r+r-) 
(bb)(bb), (bb,r+r-,_cc,gg)(r+r-), 

(r+r _ )(bb, r+r", cc,gg) 
(bb)(bb), (bb,r+r~,_cc, gg) (r+r~), 

(r+r-)(bb, r+r", cc,gg) 


189 176.2 
192 - 202 236.7 
199 - 209 217.2 


65 - 95 
60 - ^s/2 
75 - 


[32] 
[33] 
[11,34] 



Table 7: Summary of the ALEPH searches for the pair production process e + e — > Ti-zHi. The 
searches are restricted to |m^ 2 — m^J less than about 20 GeV/c 2 . 



30 





(GeV) 


£ (pb" 1 ) 


Mass ranges (GeV/c 2 ) 


Ref. 


e+e~^ Hi ( )( ) 








r u JJl c c j A 4 /- 1 )i v v A dll J/ 


91 


2.5 


< 0.21 


[361 




91 


0.5 


U.zl 


— z 


[371 


( iet^ fe^e - ^ 


91 


0.5 


1 

1 — 


on 


[371 

[O 1 j 


(ipf iptVl?+/ ,_ 7/7/ 


91 


3.6 


1Z — 


^n 




yeu jeLj^e &,[/,[/,, vv ) 


91 


33.4 






[391 


(hb)(anv) (Y+T^fnnl 


161,172 


19.9 


40- 


- 80 


[401 




183 


52.0 


45- 


- 95 


[411 


fbbVanvl (T^~T~\(nn\ 


189 


158.0 


65 - 


100 


[421 


(bb)(any) 


192-209 


452.4 


12 - 


120 


[43,44] 


(r+T - ) (qq) 


192-209 


452.4 


45- 


120 


[43,44] 


(qq, gg)(qq, ^ , e + e _ , ,u + aO 


189-209 


610.4 


4- 


116 


[45] 


e+e"^ H 2 Z^ (HiHi)Z^ (...)(...) 




mn 2 


mm 




(any)(qq) 


91 


16.2 


12-70 


< 0.21 


[46] 


(V°V°)(any but r+r") 


91 


9.7 


0.5-55 


< 0.21 


[46] 


(77) (any) 


91 


12.5 


0.5 - 60 


< 0.21 


[46] 


(4 prongs) (any) 


91 


12.9 


0.5 - 60 


0.21 - 10 


[46] 


(hadrons)(z^) 


91 


15.1 


1 - 60 


0.21 - 30 


[46] 


(t + t~t + t~)(z/;7) 


91 


15.1 


9-73 


3.5-12 


[46] 


(any)(qq, vv) 


161,172 


20.0 


40-70 


20-35 


[40] 


(bbbb)(qq) 


183 


54.0 


45 - 85 


12-40 


[41] 


(bbbb, bbcc, cccc)(qq) 


192-208 


452.4 


30 - 105 


12 - 50 


[43,44] 


(cccc) (qq) 


192-208 


452.4 


10 - 105 


4-12 


[47] 



Table 8: List of the DELPHI searches for the Higgsstrahlung processes e + e — > HiZ and 7Y 2 Z. 
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(GeV) 


C (pb- 1 ) 


Mass ran£ 


?es (GeV/c 2 ) 


Ref. 


e+e"^ n 2 n 1 ^ (...)(...) 










4 prongs 


91 


5.3 


0.2 - 10 


0.2 - 10 


[39] 


(t + t~ )(hadrons) 


91 


0.5 


4-35 


4-35 


[48] 


(r+r-) (jet jet) 


91 


3.6 


25-42 


25-42 


[49] 


(bb)(bb), (bb)(cc) 


91 


33.4 


15-46 


15-46 


[38] 


r + T~bb 


91 


79.4 


4-70 


4-70 


[47] 


bbbb 


91 


79.4 


12-40 


20 - 70 


[50] 


bbbb 


133 


6.0 


40-68 


35-73 


[51] 


bbbb, r+T'bb 


161,172 


20.0 


40-70 


35-75 


[40] 


bbbb, r+r-bb 


183 


54.0 


50 - 80 


25 - 105 


[41] 


bbbb, r+r-bb 


189 


158.0 


65 - 90 


40 - 115 


[42] 


T + r~bb 


192-208 


452.4 


50 - 100 


60 - 150 


[43,44] 


bbbb 


192-208 


452.4 


12 - 100 


40 - 190 


[43,44] 


T + r - r + r - 


189-208 


570.9 


4-90 


4- 170 


[50] 


bbbb 


189-208 


610.2 


12 - 70 


30 - 170 


[50] 


quarks or gluons 


189-208 


610.4 


4-170 


4-170 


[45] 


e+e"^ 7i 2 'H 1 ^ (...)(•••) 




m n 2 


mm 




(77) (77) 


91 


12.5 


0.5 - 60 


< 0.21 


[46] 


(4 prongs) (2 prongs) 


91 


12.9 


0.5 - 60 


0.21 - 10 


[46] 


(hadrons) (hadrons) 


91 


15.1 


1-60 


0.21 - 30 


[46] 




91 


15.1 


9-60 


3.5 - 12 


[46] 


(any) (any) 


161,172 


20.0 


40-70 


20-35 


[40] 



Table 9: List of the DELPHI searches for the pair production process e + e — > H2H1. 
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(GeV) 


C (pb- 1 ) 


Mass ranges (GeV/c 2 ) 


Ref. 


e+e-^HiZ^ (...)(...) 








(bb)(any),(r + r~)(qq) 


189 


176.4 


60 - 100 


[52] 


(bb)(any),(r+r-)(qq) 


192 - 202 


233.2 


60 - 110 


[53] 


(bb)(any),(r + r _ )(qq) 


203 - 209 


217.3 


60 - 120 


[54] 


(bb, cc, gg)(any) 


189 


176.4 


60 - 100 


[55] 


(bb, cc, gg)(any) 


192 - 202 


233.2 


60 - 110 


[55] 


(bb, cc, gg)(any) 


204 - 209 


214.5 


60 - 120 


[55] 


e+e-^ H 2 Z^ (WiWi)Z-v (...)(...) 




~mn 2 m Hi 




bb,cc,gg)(qq) 


189 - 209 


626.9 


30 - 85 10 - 42 


[56] 



Table 10: List of the L3 searches for the Higgsstrahlung processes e + e — > HiZ and H.{L. 





(GeV) C (pb- 1 ) 


Mass ranges (GeV/c 2 ) 


Ref. 


e+e"^ n 2 n 1 ^ (...)(...) 








(bb)(bb), (bb)(r+r-), (r+r") (bb) 
(bb)(bb), (bb)(r + r-), (r+r") (bb) 
(bb)(bb), (bb)(r+r-), (r+r") (bb) 


189 176.4 
192 - 202 233.2 
204 - 209 216.6 


50 - 95 50 - 95 
50 - 105 50 - 105 
50 - 110 50 - 110 


[57] 
[58] 
[56] 



Table 11: List of the L3 searches for the pair production process e + e — > T-i^Hx- 
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^(GeV) 


£(pb^) 


Mass ranges (GeV/c 2 ) 




rtl£—> (...) (...) 




mm 




(bb)(qq) 


lbl-172 


on a 

zU.4 




Ad on 

40 — 80 




rcn ftp.] 
[oy, duj 


(bbj(qqj 


loo 


54.1 




/in nr 

40 — 95 




[blj 


(bb)(qq) 


189 


1 TO 1 

l/Z.l 




/in 1 nn 

40 — 100 




[b2\ 


(bb)(qq) 


1 no orifi 


A 1 O 

4zl.z 




on 1 on 

80 — 120 




[06\ 


(00) (W J 


lol— 17z 


on /i 

zU.4 




cr n Tn 

o(J — 71) 




[ou, DUJ 


(00) (W ) 


loo 


£0 n 

00. y 




r n of; 

ou — yo 




[blj 


(OOJ (1/1/ J 


1 en 


171 /I 
1 ( 1.4 




en 1 nn 

OU — 1UU 




[bZJ 


(00) (Vis J 


1 no onn 

iyz— zuy 


a 1 n n 




on 1 on 

oU — IzU 




[bc5J 


(00)(T T ), (T T jl^qqj 


iDl— 1 / Z 


on /i 
ZU.4 




on ciK 

ou — yo 




r^o fini 

[(Jc/j UUJ 


(00)(T T ), (T T J^qq; 


loo 


00. / 




on 1 nn 
oU — 1UU 




[fill 

[01 j 


(00)(T T ), (T T J^qq; 


1 on 
loy 


loo. / 




on 1 nn 
oU — 1UU 




[DZJ 


(00)(T T ), (T T jl^qqj 


1 no onn 

iyz— zuy 


A~l 7 A 

41 ( .4 




on 1 on 
oU — IzU 




[OOJ 


(oo)(e e ), (OD)(fi \i ) 


loo 


n n 

00. y 




/?n 1 nn 
oU — 1UU 




[Olj 


(00)(e e ), (00)([i ) 


1 on 
lo9 


1 7n n 
1 /U.U 




Tn 1 nn 
/U — 1UU 




[b/j 


(00)(e e ), (OD)([i ) 


1 no onn 
lyz— zuy 


/HQ Q 
4l0.0 




/in 1 on 
4U — IzU 




[b^j 


(qq\ gg)( r+r_ ) vv), (t + t~ )(qq) 


91 


4b. 3 




n Tn 

— 7U 




[04, OOJ 


(qq> gg)(e + e", ^ + ^~) 


91 


4b. 3 




on Tn 

20 — 7U 




[04, OOJ 


(&ny)(e e , /x fi ) 


iDl— 1 t I 


on yj 

ZU.4 




or oa 

00 — 8U 




[oy, ouj 


(qq ; gg)(qq) 


189 


174.1 




60 - 100 




[OOJ 


iqq> ggjiqqj 


192-209 


424.2 




60 - 120 




[°'J 


inn <T<T 1 / 7/T/l 

iqq> ggjv^^j 


189 


171.8 




30 - 100 




[OOJ 


I nn ffrrll 7/7/I 

iqq> ggjv^^j 


192-209 


414.5 




30-110 




L° * J 


iqq> eg a* T )■> \ T T Jiqqj 


189 


168.7 




30 - 100 




[OOJ 


(hhi &&)( t T )■> \ T T Jiqqj 


192-209 


418.9 




60-115 




[Of J 


iqq> b&)(v v , n n ) 


189 


170.0 




70 - 100 




[OOJ 


iqq> ggjv e e > m r ) 


192-209 


422.0 




60 - 120 






o+p-^ 7-/„7^ (7-1 ,7-1 A7 ^ ( \( \ 
c e * n.2t-> * (rL\r\\)Li > (...)(...) 




mn 2 mm 




(qqqq)(^i/) 


91 


46.3 


10- 


75 


-35 


[04, OOJ 


i DDDD Jiqqj 


183 


54.1 


40- 


80 10.5 - 38 


[blj 


iDDDbjiqqj 


189 


172.1 


40- 


100 10.5 - 48 


te9l 
[02J 


l DDDD Jiqqj 


192-209 


421.2 


80- 


120 12 - 


m n j2 


n m 

[1UJ 




183 


53.9 


50 - 


95 10.5 - 


- mm/ 2 


tell 

[Olj 


(qqqq)(^) 


189 


171.4 


50 - 


100 10.5 - 


- m H , 2 /2 


[62] 


(bbbb)(W) 


199-209 


207.2 


100- 


110 12- 


m n , 2 /2 


[10] 


(bbbb) (r+r") 


183 


53.7 


30- 


100 10.5 - 


- m H2 /2 


[61] 


(bbbb) (t + t ) 


189 


168.7 


30- 


100 10.5 - 


- ™W 2 


[62] 


(bbbb, bbr + T~, t + t~t + t~) 














(uu, e + e~, (i + fJ>~) 


189-209 


598.5 


45- 


90 2- 


- 10.5 


[68] 



Table 12: List of the OPAL searches for the Higgsstrahlung processes e + e — > 7YiZ and 



34 





V^(GeV) 


£(pb-') 




TVTass ranpfps 


(GeV/c 2 ) 


Ref. 


n 2 n 1 ^ (...) (...) 










Ref. 


\pu )\DD ) 




loU— loD 


^ 9 
O.Z 


V 


— so — i 


A — — ^0 

i—\ — U «JU 


[601 


(DU ) (DD ) 




i ri 

lul 


inn 


y 


- sn _ i ^n 


A — — 60 

— i — u uu 


[59 601 

Ui/ 5 UU 


(DD )(DD ) 




1 79 
1 / Z 


i n /i 

1U.4 


y; 


— sn _ 1 so 

OU -LOU 


A — — 60 

t—A — U UU 


T5Q 601 

Uf ; UU 






1 SQ 
loo 


1 

04.1 


y; 


— sn _ i ^n 


A — — 60 

— i — u uu 


[611 
L U1 J 


(DD )(DD ) 




1 CO 

loy 


1 79 1 
1 ( Z.l 


y 


— 80 — i 

— OU lOU 


A — — 70 

i-4 — U / U 


[691 

[DZj 


\V)U ) (DD ) 




1 Q9 

iyz 


98 Q 

zo.y 


y 


— OO ±00 


A — — 70 

t—A — U 1 U 


[101 


(DD )(DD ) 




iyo 


/ 4.8 


y 


— sn — i S7 

— OU -LO 1 


A — n — 70 


[101 

[1UJ 


(DD )(DD ) 




zuu 


77 9 
/ I .Z 


y 


— 8f) _ 1 Q1 

OU ± t? -L 


A — — 70 

£—* — U 1 U 


[101 






zuz 


OD. 1 


y 


— OU J-c/O 


A — — 70 

i-i — u / U 


[1 nl 

[1UJ 






1 QQ_9nQ 
lyy zuy 


9D7 S 


y 


— 1 90 _ 1 qn 


A — — 70 

l— \ — U i u 


[101 


(UU )(UU J 




1 QQ 90Q 

iyy zuy 


9n7 ^ 

ZU 1 .O 


y 


- i nn _ 1 4f) 

— 1UU ±^L\J 


A — 60 — 1 00 

i— 1 — UU 1UU 


[101 

[lUj 


\UU J \ 1 1 J j \ 1 1 


)Cbb) 
i\ uu ) 


161 


10 




40 — 160 


52 - 160 


[59 601 


(bb)(r + T~) (t + t~ 

\UU J \ 1 1 J j \ 1 1 


)<bb) 
l\ uu ) 


1 72 


1 4 




37 — 160 


28 — 160 


59 60 


(bb)(-r + T~) (t + t~ 

1 V U J 1 / / / ; \ ' ' 


)(bb) 
)(UU) 


1 83 


53 7 

U(J. 1 


y 


= 70 - 170 


A = - 70 


[611 


(bb)(r + T~) (t + t~ 


)(bb) 


1 8Q 


1 68 7 


y 


= 70 - 190 


A = - 90 


[621 

L w J 


fhhVr+r - ") (V+r - 


)(bb) 
>\ UD ) 


1Q2 


98 7 


y 


= 10 - 174 


A = - 182 


[101 

L J 


(hbVr+r - ^ fr+r - 


)<bb) 
l\ uu ) 


1 96 


74 7 


y 


= 10 - 182 


A = - 191 


[101 

L J 


(bb)(-r + T~) (t + t~ 

1 V U II / / y ; \ ' ' 




2DD 

Z/UU 


74 8 


y 


= 10 - 182 


A = - 191 


[101 


(bbVr+r _ ) ('r+r - 




202 


35.4 


y 


= 10 - 174 


A = - 182 


[101 


(bb)(r+r"), (r+r~ 


)(bb) 


199-209 


203.6 


y 


= 70 - 190 


A = - 90 


[101 

L J 


(qq)(r+r~), (r+r - ) 


(qq) 


91 


46.3 




12 - 75 


10 - 78 


[64, 65] 


e+e"^ H 2 fti^ 


















i-> (•••)(•••) 








mn 2 


m«i 




(bbbb)(bb) 




91 


27.6 




40-70 


5-35 


[64, 65] 


(bbbb)(bb) 




130-136 


5.2 




55 - 65 


> 27.5 


[60] 


(bbbbj(bb) 




161 


10.0 




55 - 65 


> 20.0 


[59,60] 


(bbbb)(bb) 




172 


10.4 




55-65 


25-35 


[59,60] 


(bbbb)(bb) 




183 


54.1 




30-80 


12 - 40 


[61] 


(bbbb)(bb) 




189 


172.1 




24-80 


12-40 


[62] 


(bbbb)(bb) 




199-209 


207.3 


s 


= 90 - 200 


A = 40 - 160 


[10] 


6r, 4r2q, 2r4q 


91 


46.3 




30-75 


4-30 


[64, 65] 



Table 13: List of the OPAL searches for the pair production process e+e — > I-L2H1 ■ The symbols 
E and A stand for the mass sum m^ 2 + and mass difference |m^ 2 — m-^J. 
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Appendix B: Limits on topological cross-sections 



The tables presented below summarise the 95% CL upper bounds, as a function of the 
Higgs boson masses, of the scaling factor S 95 defined in the text (see Eq. 16). Tables ITU IT51 
and 1161 refer to final-state topologies arising from the Higgsstrahlung processes e + e - ^ TiiZ 
and e + e~— > i^H%— > 7ii7ii)Z; Tables ITS1 to 1711 refer to those arising from the pair production 
processes e + e~-^ 'H'ftii and e + e~— > r H.{H\) r H.\. The corresponding figures, showing the 

same results, are mentioned in the table captions. 
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(GeV/r 2 ) 


fa"! 




(r) 


(GeV /r 2 ) 


W 


\ u ) 




1 9 


0904 

U.U^U^± 


n m ^4 

U.UXO^i 


n 009^ 


uu 


U.UiOU 


D D91 s 


D D987 


1 4 


U.UX 1 u 


n m 4^ 


n dsqq 

u.uoyy 


fi8 
uo 


n 09^^ 

U.UiOU 


D D91 S 


D D9R7 


1 

±u 


U.UXOO 






70 


D D971 

U.U^ / L 


n n94fi 

u.u^^u 


D D9X7 


1 8 


n m ^o 

U.UXOU 


n m ^1 


U.UjOO 


79 


D9Q1 


D D974 

U.U^ / 


D D971 

U.U^ / L 




U.UXOU 


u.uioy 


U. 1UUU 


74 


n 0^90 


n n^ni 


D D9Q7 


99 


n m 77 

U.UX ( / 


U.UXOU 


n i 080 

U. 1UOU 


7fi 


D 0491 


u.uoou 


n n^i 


94 


n m Q4 


n m 74 

U.Ul 


n 1 1 1 n 

U.111U 


78 


o n4fiQ 


n 0494 


u.uoou 


9fi 


D 0907 
u.u^u / 


U.UXOU 


n 1 1 4n 


80 
ou 


u.u^oo 


n D41 n 

u.u^xu 


n n^i 

U.UOIU 


98 
zo 


n 099*} 


u.uiyo 


n 1 1 1 n 

U.111U 


89 


D 04R7 

U.U^±U 1 


D 047^ 


D 0981 

U.UiOl 


^0 


n 0901 


n m si 

u.uxox 


u.uoyo 


84 

0^± 


n o^q 

u.uooy 


n n^s^ 

U.U(JO(J 


D D999 


^9 
oz 


u.uxyo 


U.UX I o 


fl 07Qfi 


8fi 
ou 


D 07R9 
U.Ul u^ 


n nsi r 


D 09^7 

U.U^O 1 




n m qi 


o m 79 

U.UX I £j 


D 0R89 


88 
oo 


n 1 1 9 


niis 

U. -L XO 


D D9Qfi 

U.U^c/U 


OU 


D941 


n m 87 

u.uxo t 


U.UUOO 


QO 

yu 


n 1 ^ 

U. XOO 


n i ^9 

U. LO£ 


n n^i 

U.UOOl 


38 
oo 


D D9QQ 


n 09"^ 

u.u^oo 


U.UU04 


Q9 


n 1 7Q 


17^ 


U.UOO^i 




u.uooo 


D 09fi7 


U.UU1J 


Q4 

y^± 


D 99Q 


n 9i 4 


D D4Q1 


49 


u.uou / 


D D9Q7 




yu 


n 9^0 


D 99D 


U.UO / u 




fl 0^78 


n 0^1 n 

U.UO-LU 


n o^Q4 


Q8 

yo 


n 9^8 

U.^OU 


U.iOO 


u.uouo 


4fi 
^u 


n 0^87 


D 0^98 


n n^79 


1 00 

1UU 


n 944 


n 91 (\ 

U.Z1U 


u.uooz 


48 


0.0391 


0.0337 


0.0575 


102 


0.237 


0.216 


0.0588 


50 


0.0363 


0.0316 


0.0445 


104 


0.255 


0.227 


0.0704 


52 


0.0386 


0.0344 


0.0454 


106 


0.263 


0.223 


0.0896 


54 


0.0387 


0.0349 


0.0464 


108 


0.266 


0.227 


0.110 


56 


0.0384 


0.0360 


0.0403 


110 


0.297 


0.244 


0.144 


58 


0.0390 


0.0367 


0.0427 


112 


0.435 


0.343 


0.212 


60 


0.0398 


0.0365 


0.0456 


114 


0.824 


0.640 


0.410 


62 


0.0293 


0.0264 


0.0444 


116 


1.41 


1.79 


1.79 


64 


0.0278 


0.0258 


0.0394 











Table 14: The 95% CL upper bound, S 95 , obtained for the normalised cross-section (see text) of 
the Higgsstrahlung process e + e - ^ 7~CiZ, as a function of the Higgs boson mass. The numbers 
listed in this table correspond to the observed limit (full line) in Figured which is reproduced 
from Ref. [3]. In the columns labelled (a) the Higgs boson is assumed to decay as in the 
Standard Model; in columns (b) it is assumed to decay exclusively to bb and in columns (c) 
exclusively to t + t~ . 
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rn-17 (GeV/r 2 ) 


m ni (GeV/c 2 ) 




1 n 


1 ^ 

10 


90 

zu 


9"i 

ZO 


oU 


00 


4U 


4:0 


OU 


OO 


90 

zu 


U.UZU 




















9^ 
zo 


U.UZO 






















U.UO I 


U.U40 


















OO 


0.048 


0.042 


















AD 


0.053 


0.056 


0.051 
















^±(j 


0.066 


0.059 


0.046 
















50 


0.087 


0.058 


0.048 


0.049 














55 


0.11 


0.055 


0.050 


0.050 














60 


0.29 


0.103 


0.094 


0.094 


0.053 












65 


0.30 


0.099 


0.091 


0.088 


0.084 












70 


0.25 


0.098 


0.097 


0.095 


0.083 


0.059 










75 


0.34 


0.11 


0.10 


0.11 


0.10 


0.096 










80 


0.39 


0.13 


0.14 


0.14 


0.13 


0.12 


0.13 








85 


0.52 


0.20 


0.20 


0.20 


0.21 


0.19 


0.18 








90 


> 1 


0.23 


0.23 


0.23 


0.27 


0.26 


0.24 


0.28 






95 


> 1 


0.29 


0.27 


0.29 


0.31 


0.29 


0.28 


0.30 






100 


> 1 


0.30 


0.29 


0.31 


0.30 


0.27 


0.28 


0.29 


0.29 




105 


> 1 


0.27 


0.32 


0.36 


0.40 


0.36 


0.31 


0.35 


0.35 




110 


> 1 


0.44 


0.54 


0.55 


0.96 


0.97 


> 1 


> 1 


0.89 


> 1 



Table 15: The 95% CL upper bound, S95, obtained for the normalised cross-section (see text) 
of the Higgsstrahlung cascade process e + e~— > (H.2—> r H.{H\fL^> (bbbbjZ, as a function of the 
Higgs boson masses m-^ 1 and m^. The numbers correspond to the contours shown in Figure^ 

to- 
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15 


0.033 
















20 


0.048 


0.32 














25 


0.070 


0.076 














30 


0.10 


0.11 


0.38 












35 


0.18 


0.19 


0.51 












40 


0.22 


0.22 


0.40 


0.39 










45 


0.30 


0.31 


0.49 


0.49 










50 


0.18 


0.38 


0.66 


0.66 


0.63 








55 


0.18 


0.37 


0.68 


0.69 


0.68 








60 


0.20 


0.38 


0.95 


0.96 


0.96 


0.94 






65 


0.20 


0.38 


> 1 


> 1 


> 1 


> 1 






70 


0.21 


0.43 


> 1 


> 1 


> 1 


> 1 


> 1 




75 


0.19 


0.46 


> 1 


> 1 


> 1 


> 1 


> 1 




80 


0.20 


0.44 


0.83 


0.83 


0.83 


0.83 


0.84 


0.84 


85 


0.25 


0.56 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 



Table 16: The 95% CL upper bound, S95, obtained for the normalised cross-section (see text) 
of the Higgsstrahlung cascade process e + e~— > (H.2-* HiHi)Z^ (t + t~t + t~ )Z, as a function 
of the Higgs boson masses and m-ft 2 . The numbers correspond to the contours shown in 
Figured (b). 
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(GeV/c 2 ) 


(a) 


(b) 


(c) 


(d) 


(GeV/c 2 ) 


(a) 


(b) 


(c) 


(d) 





0.0237 


0.0237 


0.0237 


0.0237 


105 


0.0243 


0.0213 


0.0354 


0.0300 


5 


0.0238 


0.0238 


0.0238 


0.0238 


110 


0.0297 


0.0250 


0.0418 


0.0313 


10 


0.0242 


0.0242 


0.0242 


0.0242 


115 


0.0472 


0.0387 


0.0484 


0.0332 


15 


0.0248 


0.0248 


0.0248 


0.0248 


120 


0.0682 


0.0599 


0.0409 


0.0348 


20 


0.0255 


0.0255 


0.0255 


0.0255 


125 


0.0676 


0.0542 


0.0493 


0.0387 


25 


0.0266 


0.0266 


0.0266 


0.0042 


130 


0.0688 


0.0541 


0.0524 


0.0429 


30 


0.0054 


0.0054 


0.0018 


0.0043 


135 


0.0618 


0.0478 


0.0571 


0.0604 


35 


0.0044 


0.0041 


0.0018 


0.0043 


140 


0.0669 


0.0524 


0.0660 


0.0665 


40 


0.0029 


0.0026 


0.0021 


0.0048 


145 


0.0600 


0.0540 


0.0506 


0.0739 


45 


0.0033 


0.0030 


0.0021 


0.0051 


150 


0.0798 


0.0726 


0.0591 


0.0847 


50 


0.0036 


0.0034 


0.0017 


0.0055 


155 


0.0967 


0.0895 


0.0696 


0.0995 


55 


0.0043 


0.0042 


0.0016 


0.0067 


160 


0.136 


0.125 


0.0847 


0.118 


60 


0.0055 


0.0057 


0.0016 


0.0083 


165 


0.179 


0.122 


0.175 


0.144 


65 


0.0073 


0.0070 


0.0010 


0.0097 


170 


0.323 


0.237 


0.234 


0.188 


70 


0.0097 


0.0106 


0.0021 


0.0117 


175 


0.352 


0.294 


0.245 


0.269 


75 


0.0142 


0.0163 


0.0029 


0.0134 


180 


0.765 


0.596 


0.408 


0.391 


80 


0.0203 


0.0227 


0.0043 


0.0165 


185 


0.838 


0.702 


0.582 


0.700 


85 


0.0357 


0.0383 


0.0101 


0.0198 


190 


1.04 


0.855 


0.764 


1.07 


90 


0.0527 


0.0522 


0.0292 


0.0247 


195 


1.93 


1.81 


1.10 


2.88 


95 


0.0520 


0.0493 


0.0400 


0.0266 


200 


6.97 


6.47 


3.49 


5.29 


100 


0.0298 


0.0257 


0.0370 


0.0283 













Table 17: The 95% CL upper bound, S 95 , obtained for the normalised cross-section (see text) 
of the pair production process e + e~-^ H.2H1, as a function of the Higgs boson mass sum 
m^i + rn-H.2 ■ Tn e bounds are given for the particular case where my_ 2 and mu-y are approxi- 
mately equal. This occurs, for example, in the CP-conserving MSSM scenario m^-max for tan f3 
greater than 10 and small m-^ 2 (= m\). The numbers listed in this table correspond to the four 
plots in Figured (see the corresponding labels). For m ni + m H2 less than 30 GeV/c 2 , the 
bounds are derived from the measured decay width of the Z boson, see Section 3.2. Columns 
labelled (a): the Higgs boson decay branching ratios correspond to the m^-max benchmark sce- 
nario with tan (3=10, giving 94% for Hi^ bb, 6% for H\-+ t + t~, 92% for H 2 ^ bb and 8% for 
t + t~ ; columns (b): both Higgs bosons are assumed to decay exclusively to bb; columns 
(c): one Higgs boson is assumed to decay exclusively to bb only and the other exclusively to 
t + t~ ; columns (d): both Higgs bosons are assumed to decay exclusively to t + t~ . 



40 



m W2 (GeV/c z ) 










(GeV/ 


(?) 










10 


15 


20 


25 


30 


35 


40 


45 


50 


15 


> 1 
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0.011 
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0.023 
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0.028 
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0.043 
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0.029 


0.026 
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0.048 


0.067 


0.041 
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> 1 
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0.026 


0.042 


0.044 


0.069 


0.043 


0.035 


0.028 


55 


> 1 
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0.050 


0.038 


0.030 


60 


> 1 


0.075 
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0.098 


0.051 


0.051 


0.050 


0.044 


0.039 


65 


> 1 


0.14 


0.13 


0.10 


0.065 


0.064 


0.070 


0.068 


0.069 


70 


> 1 


0.20 


0.16 


0.11 


0.072 


0.074 


0.066 


0.072 


0.071 


75 


> 1 


0.23 


0.13 


0.14 


0.076 
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0.083 


0.066 


0.093 


80 


> 1 


0.26 


0.19 


0.12 


0.078 


0.089 
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0.064 


0.093 


85 


> 1 


0.26 


0.17 


0.13 


0.095 


0.080 


0.070 


0.071 


0.10 


90 


> 1 


0.18 


0.13 


0.11 


0.073 


0.070 


0.076 


0.081 


0.13 


95 


> 1 


0.20 


0.13 


0.095 


0.073 


0.078 


0.081 


0.11 


0.15 


100 


> 1 


0.21 


0.12 


0.092 


0.085 


0.091 


0.12 


0.16 


0.18 


105 


> 1 


0.16 


0.13 


0.12 


0.11 


0.13 


0.18 


0.20 


0.20 


110 


0.297 


0.16 


0.14 


0.15 


0.14 


0.17 


0.20 


0.20 


0.19 


115 


0.338 


0.22 


0.20 


0.20 


0.18 


0.20 


0.21 


0.21 


0.23 


120 


0.355 


0.28 


0.25 


0.23 


0.22 


0.22 


0.23 


0.27 


0.36 


125 


0.409 


0.29 


0.26 


0.25 


0.22 


0.25 


0.29 


0.40 


0.51 


130 


0.494 


0.35 


0.32 


0.32 


0.24 


0.32 


0.46 


0.57 


0.72 


135 


0.617 


0.44 


0.42 


0.42 


0.36 


0.51 


0.67 


0.84 


0.98 


140 


0.696 


0.57 


0.53 


0.66 


0.62 


0.83 


0.97 


> 1 


> 1 


145 


0.811 


0.73 


0.80 


> 1 


0.94 


> 1 


> 1 


> 1 
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0.082 


0.078 
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0.10 
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0.098 
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0.11 


0.14 


0.14 












85 


0.12 


0.14 


0.16 


0.15 


0.21 










90 


0.17 


0.16 


0.17 


0.25 


0.24 


0.41 








95 


0.19 


0.21 


0.30 


0.35 


0.44 


0.51 


0.64 






100 


0.21 


0.31 


0.39 


0.42 


0.43 


0.74 


> 1 


> 1 




105 


0.32 


0.42 


0.55 


0.53 


0.90 


> 1 


> 1 


> 1 


> 1 


110 


0.47 


0.55 


0.63 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


115 


0.56 


0.65 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


120 


0.64 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


125 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 



Table 18: The 95% CL upper bound, S 95 , obtained for the normalised cross-section (see text) 
of the pair production process e + e~— > 'H^Hx^ bbbb, as a function of the Higgs boson masses 
rriHi an d Tn-H 2 . The numbers correspond to the contours shown in Figured (a). 
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0.0065 
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0.011 
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0.018 
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0.17 
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115 


0.14 


0.20 


0.27 


0.36 
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> 1 


> 1 
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0.19 
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0.83 


> 1 


> 1 


> 1 
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> 1 
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0.53 
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> 1 


> 1 


> 1 


> 1 
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> 1 


> 1 
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0.46 


0.85 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


135 


0.89 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 



Table 19: The 95% CL upper bound, S95, obtained for the normalised cross-section (see text) 
of the pair production process e + e~— > 7i 2 7ii^ t + t~t + t~ , as a function of the Higgs boson 
masses and m-^ 2 . The numbers correspond to the contours shown in Figure^ (b). 
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Table 20: The 95% CL upper bound, S95, obtained for the normalised cross-section (see text) of 
the pair production cascade process e + e~— > (% 2 ^ 'H\H\)'H.\-^ (bbbbjbb, as a function of the 
Higgs boson masses m Wl and m n2 . The numbers correspond to the contours shown in Figure® 
(a). 



43 



m H2 (GeV/c 2 ) 


m Wl (GeV/c 2 ) 




5 


10 


15 


20 


10 


0.0006 








15 


0.0016 








20 


0.0017 


0.011 






25 


0.0018 


0.0019 






30 


0.0021 


0.0021 


0.013 




35 


0.0024 


0.0025 


0.017 




40 


0.0009 


0.0016 


> 1 


> 1 


45 


0.0010 


0.0019 


> 1 


> 1 


50 


0.0013 


0.0023 


> 1 


> 1 


55 


0.0017 


0.0029 


> 1 


> 1 


60 


0.0024 


0.0043 


> 1 


> 1 


65 
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Table 21: The 95% CL upper bound, S 95 , obtained for the normalised cross-section (see text) of 
the pair production cascade process e + e~ — > (H<i-^ TCiTCi)TCi^> (t + t~t + t~)t + T ; 3jS Si function 
of the Higgs boson masses and m^ 2 . The numbers correspond to the contours shown in 
Figure® (b). 
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Figure 1: Contours of the observed p-values, 1 — CL&, indicating the statistical significances 
of local excesses in the data. Plots (a) and (b) refer to the CP-conserving MSSM benchmark 
scenario m^-max and plots (c) and (d) to the CP-violating scenario CPX. For each scenario, the 
parameter space is shown in two projections. Regions which are not part of the parameter space 
(labelled "Theoretically Inaccessible") are shown in light-grey or yellow. In the medium-grey 
or light-green regions the data show an excess of less than one standard deviation above the 
expected background. Similarly, in the dark-grey or dark-green regions the excess is between 
one and two standard deviations while in the darkest-grey or blue regions it is between two 
and three standard deviations. In plots (c) and (d), two small regions with excesses larger than 
three standard deviations are shown in white. The dashed lines show the expected exclusion 
limit at 95% CL. The hatched areas represent regions where the median expected value of CL S 
in the background hypothesis is larger than 0.4; apparent excesses in these regions would not 
be significant. 
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Figure 2: The 95% CL upper bounds, S95 (see text), for various topological cross-sections 
motivated by the Higgsstrahlung process e + e~^ TiiZ, as a function of the Higgs boson mass 
(the figure is reproduced from Ref. [3]). The full lines represent the observed limits. The 
dark (green) and light (yellow) shaded bands around the median expectations (dashed lines) 
correspond to the 68% and 95% probability bands. The horizontal lines correspond to the 
Standard Model cross-sections. In part (a ) the Higgs boson decay branching ratios are assumed 
to be those predicted by the Standard Model; in part (b ) the Higgs boson is assumed to decay 
exclusively to bb and in part (c) exclusively to t + t~ . 
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Figure 3: Contours of the 95% CL upper bound, S95 (see text), for various topological cross- 
sections motivated by the Higgsstrahlung cascade process e + e - ^ (H.i~^ 'H\H\JL, projected 
onto the (m n2 , m Hl ) plane. The scales for the shadings are given on the right-hand side of each 
plot. In plot (a) the Hi boson is assumed to decay exclusively to bb and in plot (b) exclusively 
to t + t~ ; in plot (c) it is assumed to decay with equal probabilities to bb and to t + t~ . 
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Figure 4: The 95% CL upper bounds, S 95 (see text), for various topological cross-sections 
motivated by the pair production process e + e~-^ H.<{Hi. The bounds are obtained for the 
particular case where m^ 2 and are approximately equal. Such is the case, for example, in 
the CP-conserving MSSM scenario m^-max for tan/? greater than 10 and small m-H 2 (= m\). 
The abscissa is the sum of the two Higgs boson masses. The full lines represent the observed 
limits. The dark (green) and light (yellow) shaded bands around the median expectations 
(dashed lines) correspond to the 68% and 95% probability bands. The curves which complete 
the exclusions at low masses are obtained using the constraint from the measured decay width 
of the Z boson, see Section 3.2. Plot (a): the Higgs boson decay branching ratios correspond 
to the m^-max benchmark scenario with tan (3=10, namely 94% Ti\-^ bb, 6% 7i\— > t + t~ , 92% 
bb and 8% H.^-^ t + t~; plot (b): both Higgs bosons are assumed to decay exclusively 
to bb; plot (c): one of the Higgs bosons is assumed to decay exclusively to bb and the other 
exclsively to t + t~ ; plot (d): both Higgs bosons are assumed to decay exclusively to t + t~ . 
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Figure 5: Contours of the 95% CL upper bound, S95 (see text), for various topological cross- 
sections motivated by the pair production process e + e~— > 'H2H1, projected onto the (m n , 2 , 
m Hi) plane. The scales in terms of the shadings are given on the right-hand side of each plot. 
In plot (a) both Higgs bosons are assumed to decay exclusively to bb and in plot (b) exclusively 
to t + t~ . In plot (c) the Ti,2 boson is assumed to decay exclusively to bb and the 7i\ boson 
exclusively to t + t~ and in plot (d) the 7i\ boson is assumed to decay exclusively to bb and the 
TC2 boson exclusively to t + t~ . The dashed lines represent the approximate kinematic limits of 
the processes. 
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Figure 6: Contours of the 95% CL upper bound, S 95 (see text), for various topological cross- 
sections motivated by the pair production cascade process e + e~— > (Hi—* HiHi)Hi, projected 
onto the (m^ 2 , ) plane. The scales in terms of the shadings are given on the right-hand 
side of each plot. In plot (a) the Hi boson is assumed to decay exclusively to bb and in plot 
(b) exclusively to t + t~ . In plot (c) the Hi boson is assumed to decay with equal probability 
to bb and to t + t~ . The dashed line in part (a) represents the approximate kinematic limit of 
the process. 
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Figure 7: Exclusions, at 95% CL (medium-grey or light-green) and the 99.7% CL (dark- 
grey or dark-green), in the case of the CP-conserving m^-max benchmark scenario, for m t = 
174.3 GeV/c 2 . The figure shows the theoretically inaccessible domains (light-grey or yellow) 
and the regions excluded by this search, in four projections of the MSSM parameters: (a): 
(m h , m^); (b): (m h , tan/3); (c): (m^, tan/3J; (d): (m R ±, tan/3). The dashed lines indicate the 
boundaries of the regions which are expected to be excluded, at 95% CL, on the basis of Monte 
Carlo simulations with no signal. In the (mh, tan (3) projection (plot (b)), the upper boundary 
of the parameter space is indicated for four values of the top quark mass; from left to right: m t 
= 169.3, 174.3, 179.3 and 183.0 GeV/c 2 . 
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Figure 8: Domains of tan (3 which are excluded at the 95% CL (light- grey or light- green) and 
the 99.7% CL (dark-grey or dark-green), in the case of the CP-conserving m^-max benchmark 
scenario, as a function of the assumed top quark mass. 
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Figure 9: Exclusions in the case of the CP-conserving m^-max benchmark scenario, variant 
(a) (see Section 2.1.1.). See the caption of Figure Q for the legend. Note the small domains 
at rrih between 60 and 75 GeV/c 2 , small mx and tan/3 < 0.9 which, although excluded at the 
95% CL, are not excluded at the 99.7% CL. 
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Figure 10: Exclusions in the case of the CP-conserving m^-max benchmark scenario, variant 
(b) (see Section 2.1.1.). See the caption of Figured for the legend. 
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Figure 11: Exclusions in the case of the CP-conserving no-mixing benchmark scenario. See the 
caption of Figure Q for the legend. Note the small domain at between 75 and 80 GeV / c 2 , 
small m>A and tan (3 < 0.7 which is not excluded at the 95% CL. 
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Figure 12: Exclusions in the case of the CP-conserving no-mixing benchmark scenario, variant 
(a) (see Section 2.1.2). See the caption of Figured for the legend. In the hatched domain 
(tan/3 < 0.7), the contributions from top and stop quark loops to the radiative corrections are 
large and uncertain. Note the small domain at between 56 and 72 GeV/c 2 , small niA and 
tan/3 < 1 which, although excluded at the 95% CL, is not excluded at the 99.7% CL. 
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Figure 13: Exclusions in the case of the CP-conserving large-fi benchmark scenario (see Section 
2.1.3). See the caption of Figure\]\for the legend. In the hatched domain (tan/3 < 0.7), the 
contributions from top and stop quark loops to the radiative corrections become large and 
uncertain; hence, no exclusions can be claimed there. 
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Figure 14: Exclusions in the case of the gluoph 
the caption of Figured for the legend. 
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bic benchmark scenario (see Section 2.1.3). See 
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Figure 15: Exclusions in the case of the CP-conserving small-a e ff benchmark scenario (see 
Section 2.1.3). See the caption of Figured for the legend. 
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Figure 16: Exclusions, at 95% CL (medium- grey or light-green) and the 99.7% CL (dark-grey 
or dark-green), for the CP-violating CPX scenario with m t = 174.3 GeV/c 2 . The figure shows 
the theoretically inaccessible domains (light-grey or yellow) and the regions excluded by the 
present search, in four projections of the MSSM parameter space: (m ni , m H2 ), (m ni , tan/3), 
(m-n 2 , tan f3) and (m H + , tan (3) . The dashed lines indicate the boundaries of the regions expected 
to be excluded, at the 95% CL, on the basis of Monte Carlo simulations with no signal. In 
each scan point, the more conservative of the two theoretical calculations, FeynHiggs or CPH, 
is used. 
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Figure 17: Cross-sections, as a function of tan j3, for some of the dominant signal processes, in 
the CP-violating scenario CPX, using the FeynHiggs calculation, with a centre-of-mass energy 
of 202 GeV, m t = 175 GeV/c 2 , and m Hl between 35 and 45 GeV/c 2 . 
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Figure 18: Exclusions, in the case of the CP-violating CPX scenario, for the two theoretical 
approaches, CPH and FeynHiggs. See the caption of Figure H^l for the legend. In part (a) the 
CPH calculation is used and in part (b) the FeynHiggs calculation. In part (c) the procedure 
is adopted where, in each scan point of the parameter space, the more conservative of the two 
calculations is used. 
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Figure 19: Exclusions, in the case of the CP-violating CPX scenario, for four top quark masses: 
m t = 169.3 GeV/c 2 , 174.3 GeV/c 2 , 179.3 GeV/c 2 and 183.0 GeV/c 2 . See the caption of 
Figure fl^l for the legend. 
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Figure 20: Exclusions, in the case of the CPX scenario with various CP-violating phases, 
arg(A) = arg(mg): 0°, 30°, 60°, 90° (the CPX value), 135° and 180°. See the caption of 
Figure [77)1 for the legend. In the hatched region in part (f) the calculations are uncertain (see 
text). 
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Figure 21: Exclusions, for the CP-violating CP X scenario with various values of the Higgs mass 
parameter /i: 500 GeV, 1000 GeV, 2000 GeV (the standard CPX value) and 4000 GeV. See 
the caption of Figure for the legend. In the hatched region in part (d ) the calculations are 
uncertain (see text). 
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Figure 22: Exclusions, for the CP-violating CPX scenario with various values of the soft SUSY- 
breaking scale M SU sy- (a): M SU sy =500 GeV (the standard CPX value); (b): M SVSY =1000 GeV 
while all other parameters are kept at their standard CPX values; (c): Msvsy=1000 GeV while 
A, mg and \x are "scaled" to 2000 GeV, 2000 GeV and 4000 GeV, respectively. See the caption 
of Figure for the legend. 
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